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Abstract 


Results  are  presented  for  the  analysis  of  methods  for 
predicting  total,  merchantable,  and  defect  volumes  of 
individual  trees.  The  methodology  implemented  in  comparing 
methods  included  regression  analysis,  analysis  of  variance, 
and  multiple  comparison  analysis. 

The  study  area  was  confined  to  four  regions  of  Alberta 
representing  a  variety  of  forested  ecosystems.  For  each  of 
the  three  types  of  tree  volume,  prediction  methods  were 
chosen  for  testing  based  on  prior  accuracy  and  popularity. 
The  methods  presently  in  use  by  the  Alberta  Forest  Service 
were  included  in  the  study. 

For  total  volume,  standard  and  local  volume  functions 
were  compared.  The  standard  volume  functions  were 
significantly  more  accurate.  Species  or  regional  differences 
were  generally  not  significantly  different,  once  differences 
due  to  changing  heights  and  diameters  were  controlled. 

For  the  present  merchantable  volume  requirements  of  the 
AFS,  both  functions  which  predict  a  merchantable  ratio  and 
simple  taper  functions  are  very  accurate.  Taper  functions 
are  slightly  more  precise. 

The  estimation  of  defect  volumes  was  assessed  using 
regression  analysis  with  discrete  variables  such  as  species, 
and  continuous  variables  such  as  age.  The  regression  model 
developed  in  this  analysis  was  significant  in  reducing  the 
variance  of  defect  percentages,  however  R2  values  were  not 
over  40  percent. 
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1 .  Introduction 


The  accurate  estimation  of  stand  and  forest  tree 
volumes  is  largely  dependent  on  the  accuracy  of  individual 
tree  volume  estimation.  In  Alberta,  the  change  from  imperial 
to  metric  units  of  measure,  coupled  with  the  implementation 
of  a  new  province-wide  forest  inventory  resulted  in  a  need 
for  new  individual  tree  volume  functions.  Also,  the  recent 
development  of  computer  technology  meant  that  more  precise 
equations  could  be  calculated. 

This  research  project  was  initiated  to  review  and 
assess  the  methods  for  estimating  total,  merchantable,  and 
defect  volumes.  The  project  was  established  with  the 
cooperation  and  support  of  the  Alberta  Forest  Service  (AFS), 
Timber  Management  Branch.  Data  were  collected  as  part  of  the 
current  Phase  III  inventory  program,  under  the  jurisdictions 
of  the  AFS  and  the  Resource  Evaluation  and  Planning  Branch 
(REAP).  Four  regions  and  major  commercial  species  of  those 
regions  were  included  in  the  study.  The  regions  represented 
a  variety  of  forested  ecosystems  across  Alberta. 

For  total  tree  volume,  or  volume  of  the  main  stem  from 
ground  to  tree  tip,  various  local  and  standard  .  volume 
functions  were  tested  for  accuracy.  In  addition,  the 
variation  between  species,  and  between  Volume  Sampling 
Regions  (VSRs)  was  assessed. 

Merchantable  volume  was  defined  as  the  volume  of  the 
main  stem  within  utilization  specifications.  Models  tested 
for  estimating  this  type  of  volume  were  either  1)  functions 
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which  estimate  the  ratio  of  merchantable  volume  over  total 
volume  or  2)  mathematical  integrations  of  functions  which 
express  stem  form. 

Two  product  types,  pulp  and  sawlog,  were  differentiated 
for  internal-  defect  volume  analysis.  A  cull  study  was  used 
to  determine  the  amount  of  saw  and  pulp  defect  in  the  sample 
trees.  The  estimation  of  each  type  of  defect  volume  was 
reviewed  by  developing  and  testing  a  regression  model,  using 
both  discrete  and  continuous  variables.  To  develop  the 
model,  the  discrete  variables,  species,  cull  suspect  class, 
and  region,  were  individually  tested  to  determine  whether 
they  significantly  improved  the  regression.  Age  and  DBH , 
continuous  variables,  were  also  tested  for  significance. 

The  assessment  of  alternative  methods  for  estimating 
each  type  of  tree  volume  was  based  on  statistical  procedures 
including  regression  analysis,  analysis  of  variance,  and 
multiple  comparison  analysis.  Methods  presently  in  use  by 
the  AFS  for  estimating  the  three  types  of  volume  were 
included  in  this  study. 


2.  Tree  Volume  Data — Source  and  Scope 


2.1  Collection  of  Tree  Volume  Data 

Tree  volumes  were  obtained  by  cutting  felled  sample 
trees  into  short  lengths  and  determining  the  volume  of  each 
section  using  geometric  formulae.  Approximately  4500  samples 
were  collected  by  the  Alberta  Forest  Service  from  four  VSRs 
of  the  province: 

VSR  3.  Rocky  Highland.  Rocky-Clearwater  Region, 

high  elevation.  662  trees. 

VSR  2.  Rocky  Lowland.  Rocky-Clearwater  Region, 

low  elevation.  1087  trees. 

VSR  4.  Centre.  Whi tecour t-Slave  Lake  Region. 

1023  trees. 

VSR  8.  Northeast.  Lac  La  Biche-Athabasca  Region. 

1795  trees. 

The  boundaries  of  the  four  VSRs  were  predetermined  by  the 
AFS  using  elevation,  climate,  and  topography  as  well  as 
other  criteria  for  differentiation  (Figure  1). 

Sample  trees  were  located  on  either  fixed  or  variable 
area  plots.  The  allocation  of  these  plots  was  systematic 
(fixed  area  plots)  or  stratified  by  AFS  timber  type 
(variable  area  plots)  (Alberta  Forest  Service,  1981b).  The 
fixed  area  plots  correspond  to  Large  Scale  Photo  plots;  the 
variable  area  plots  are  also  temporary  sample  plots  for 
stand  volume  estimation  (Alberta  Forest  Service,  1981a). 
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The  distribution  of  plots  was  assumed  to  be  random,  although 
some  bias  may  have  occurred  through  greater  sampling  of  more 
accessible  areas. 

All  trees  1.1  centimetres1  and  greater  at  DBH,  within 
the  tree  section  plot  were  felled.  Most  of  the  tree  section 
plots  were  variable  area  plots;  hence,  the  distribution  of 
diameters  is  slightly  skewed,  because  larger  diameters  have 
a  higher  probability  of  being  sampled  (Smith,  1973).  The 
advantage  to  this  type  of  selection  is  that  more  valuable 
trees  (i.e.  larger  DBH)  will  be  better  represented  in  the 
samples . 

For  all  trees  in  a  tree  section  sample  plot,  the  DBH, 
species,  and  location  parameters  were  recorded  (Alberta 
Forest  Service,  1981a).  Each  tree  was  then  felled,  and  cut 
into  short  lengths  (Figure  2).  The  first  section  was  at  0.3 
meters  above  ground  (stump  height),  the  second  section  was 
at  1.3  metres  above  ground  (breast  height),  and  the  third 
section  was  at  2.8  metres  above  ground.  Successive  cuts  were 
at  2.5  metre  intervals  up  to  the  tree  tip.  All  major  limbs 
were  cut  in  forked  trees.  For  each  section,  diameter  inside 
bark  at  the  top  of  the  section,  the  section  length,  and  the 
dimensions  of  internal  defect  were  recorded. 


’Because  both  metric  and  imperial  measures  were  taken,  the 
following  soft  conversions  were  used;  1.1  cent imet res= 
0.4inches;  1.3  metres  (breast  height)=4.5  feet;  2.5  metres 
(section  length)=8.0  feet;0.30  metres  (stump  height)=1.0 
foot . 
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4  2.5  M 


BREAST 

HEIGHT 
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(33> 
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1.5  M 


<2.5M 
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height 


<2.5  M 


2.5  M 


42.5  M 


Figure  2.  Pictorial  Representation  of  Tree  Sectioning 
Procedure 

Sections  with  internal  defect  were  further  cut  to  determine 
the  length  of  the  defect. 

Because  of  the  long  sampling  period  (from  1973  to 
1981),  two  types  of  tally  sheets  were  used  in  recording  tree 
section  data.  For  trees  sectioned  after  1977,  additional 
information  such  as  diameter  outside  bark  at  the  top  of  each 
section  was  recorded  (Appendix  I);  also,  the  advent  of 
metr i f icat ion  in  Canada  resulted  in  a  change  in  field 
records  from  imperial  to  metric  units.  Hence,  all  data  had 
to  be  standardized  into  a  similar  format  and  all  imperial 
data  converted  to  metric.  Computer  programs  were  designed  to 
reorganize  the  data,  converting  information  to  metric  where 
necessary.  (Appendix  II,  Flow  chart  of  Computer  Programs). 
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2.2  Calculation  of  Individual  Tree  Volumes 

Once  the  data  were  keypunched,  edited,  and  rewritten 
into  a  standard  metric  format,  a  FORTRAN  (Cress  et  al . , 
1970)  computer  program,  TMA  (Tree  Measurement  Analysis)  was 
developed  to  provide  the  individual  tree  volumes.  Geometric 
formulae  were  used  to  calculate  total,  merchantable,  and  net 
volumes . 


2.2.1  Total  Tree  Volume 

Total  volume  for  each  tree  was  calculated  by  summing 
section  volumes  using  Smalian's  formula  (paraboloid 
frustrum)  for  the  main  bole  sections,  conical  formula  for 
the  top  section  and  cylindrical  formula  for  the  stump 
section . 


where 


Smalian's  h  (A  +  A  )  =  V 

2  b  u 

Conical  h  (A  )  =  V 

3  b 

Cylindrical  h  (A  )  =  V 

u 


V 

i  s 

the 

volume 

of  the  section 

(m3 ) 

A 

i  s 

the 

cross-sectional 

area 

at  the 

base 

b 

of 

the 

section 

(m2  ) 

A 

i  s 

the 

cross-sectional 

area 

at  the 

top 

u 

of 

the 

section 

(m2  ) 

h 

i  s 

the 

section 

length 

(m) 

(Husch 

et  a  7 

.,  1972) 

The  use  of  the  three  formulae  was  requested  by  the  AFS .  Many 
other  studies  have  used  the  same  equations  in  calculating 
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total  tree  volume  (Burkhart,  1977;  Beagle,  1974;  Brown, 
1934);  Husch  and  others  (1972)  show  this  to  be  an  adequate 
representation  of  the  volume  of  the  tree.  Some 
overestimation  can  occur  when  Smalian's  equation  is  used  to 
calculate  volume  over  a  long  section  length.  Husch  and 
others  (1972)  cite  a  study  which  showed  a  9.0  percent 
overestimation  in  volume  using  Smalian's  equation  on  8 
(approximately  2.5  metre)  to  16  foot  lengths.  The  advantage 
in  using  Smalian's  equation,  though,  is  that  the  parameters 
required  in  calculation  are  easily  measured. 

2.2.2  Merchantable  Volume 

Merchantable  volume  is  the  volume  of  the  tree  within 
specified  utilization  limits.  Four  specifications  were 
provided  by  the  AFS  for  merchantable  volume. 


1 . 

From 

a 

0.30 

metre 

stump 

to 

a 

7  ■ 

cm  top  diameter ( i . b . ) 

2. 

From 

a 

0.30 

metre 

stump 

to 

a 

10 

cm  top  diameter ( i . b . ) 

3. 

From 

a 

0.30 

metre 

stump 

to 

a 

13 

cm  top  diameter ( i . b . ) 

4. 

From 

a 

0.30 

metre 

stump 

to 

a 

15 

cm  top  diameter ( i . b . ) 

The  TMA  computer  program  simulates  cutting  of  the  tree  at 
the  stump  (base  cutoff  point)  and  at  the  top  (top  cutoff 
point).  The  cutoff  points  specified  by  the  AFS  were  a  height 
cutoff  at  the  base  and  a  minimum  diameter  cutoff  at  the  top 
of  the  stem. 2 


2  The  TMA  program  was  designed  to  use  either  a  height  or 
diameter  at  both  the  top  and  base  cutoff  points  for  future 
use . 
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In  cases  where  the  cutoff  point  fell  within  a  section, 
interpolation  was  used  to  estimate  length  and  diameter  for 
that  sub-section.  For  the  base  cutoff  point  the  following 
formulae  were  used: 


TSEC  =  SUM  SECL  -  Base  cutoff  point(Stump  height) 

(1) 


DBASE  =  SQRT  DIB2  - 

top 


(DIB2  -  DIB2  ) 
top  base 


SECL/TSEC 

(2) 


where  DBASE  is  the  diameter ( i . b . )  at  stump  height(cm) 

DIB  is  the  d.i.b.  at  the  top  of  the  section(cm) 
top 

DIB  is  the  d.i.b.  at  the  base  of  the  section(cm) 
base 

TSEC  is  the  length  of  the  interpolated  section(m) 
SECL  is  the  section  length(m) 

SUM  SECL  is  the  sum  of  the  sections  to  the  top 
of  the  section  to  be  interpolated (m) 

Stump  height  is  in  metres 
SQRT  is  the  square  root. 


SUM  SECL 


Figure  3.  An  Illustration  of  Parameters  Used  in  Calculating 
the  Volume  Above  the  Base  Utilization  Standard 
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To  interpolate  length  to  the  top  cutoff  point,  the  following 
equation  was  used: 


(DIB2  -  DTOP 2 ) 
base 

TSEC  =  SECL  X  - - 

(DIB2  -  DIB2  ) 
base  top 


(3) 


where  DTOP  is  the  specified  minimum  top  d.i.b.(cm) 
All  other  variables  as  in  Equation  2. 


►  SECL 


TSEC 

> 


Figure  4.  An  Illustration  of  Parameters  Used  in  Calculating 
the  Volume  Below  the  Top  Utilization  Limit 


Section  volumes  below  the  stump  height,  and  above  the 


minimum  diameter  were  set  to  zero. 


2.2.3  Defect  Volume 


Internal  defect  volume  was  calculated  by  section  for 
two  product  types,  pulp  and  sawlog.  For  both  types,  the 
Alberta  Scaling  Manual  (1979)  was  used  as  a  guide  in 
determining  the  defect  volume  for  each  type.  Because  the 
pattern  of  defect  varies  greatly,  a  hand-held  electronic 
calculator  was  used  for  more  subjective  calculations. 

Internal  defect  for  sawlog  products  (sawlog  defect)  was 
calculated  by  taking  a  rectangular  cross-sect ional  area  of 
defect  in  square  metres,  and  multiplying  by  the  section 
length  (metres).  For  the  stump  section,  the  cross-sectional 
area  of  defect  at  the  top  of  the  stump  was  used  in 
calculating  the  section  defect  volume.  For  the  second  and 
third  sections,  the  larger  of  the  defect  area  at  the  base  of 
the  second  section,  or  at  the  top  of  the  third  section,  was 
projected  through  both  of  the  lengths  to  simulate  a  2.5 
metre  log.  Defect  volumes  in  subsequent  sections  were 
assessed  by  taking  the  larger  of  the  squared  defect  area  at 
the  base  or  at  the  top  of  the  section,  and  multiplying  by 
the  length  of  the  defect.  A  two  centimetre  (or  one  inch) 
allowance  was  added  to  the  actual  dimensions  of  defect 
tallied  on  the  field  records  as  allowable  loss  in  sawlog 
production  (Alberta  Forest  Service,  1979). 

For  pulp  products,  allowable  loss  was  calculated  using 


Smalian's  equation  when  defect  extended  through  the  section, 
otherwise  the  actual  length  of  the  defect  was  used  with  an 
equation  for  calculating  the  volume  of  a  paraboloid.  Again, 
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this  method  is  similar  to  log  scaling  procedures. 

For  each  of  the  defect  categories,  pulp  and  sawlog,  net 
volumes  were  calculated: 

1.  From  ground  to  tree  tip  (net  total  pulp  and  sawlog 
volumes ) . 

2.  From  0.3  metres  above  ground  to  7.0,  10.0,  13.0,  and 
15.0  cm  top  diameter  inside  bark  (net  merchantable  pulp 
and  sawlog  volumes). 

When  the  base  or  top  cutoff  point  fell  within  a  section, 
defect  volume  was  estimated  using  the  following  formulae: 

(DIB2  +  DBASE2  or  DTOP2  +  DIB2  )  X  TSEC 
top  base 

drat  =  - - - 

(DIB2  +  DIB2  )  X  SECL 


top  base 


(4) 


DEFECT  =  DRAT  X  DEFECT  (pulp/saw) 
R 


where : 

DEFECT  is  the  defect  volume  for  the  reduced  section(m3) 

R 

DRAT  is  a  ratio  of  volume  for  the  reduced  section, 
to  total  section  volume 

All  other  variables  as  defined  for  Equation  2. 

Section  net  volumes  were  calculated  by  subtracting  defect 
from  gross  volumes.  Because  the  squared  defect  method  was 
used  for  sawlog  defect,  some  sections  had  more  defect  volume 
than  gross  volume;  net  saw  volume  was  set  to  zero  for  these 
sections.  Section  net  volumes  were  totaled  to  yield  net 
total  and  merchantable  tree  volumes. 
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2.3  Distribution  of  Data 

Table  1  shows  the  distribution  of  all  samples  by 
species  and  VSR.  Only  those  VSR/species  groups  of  25  or  more 

samples  were  included  in  testing.  Even  though,  for  some 

* 

species,  the  combination  of  all  VSRs  would  have  resulted  in 
over  25  samples,  the  analysis  was  not  attempted  for  the 
following  reasons: 

1.  As  more  areas  and/or  more  species  are  combined  in  a 
sample  unit,  the  variance  increases  so  a  greater  number 
of  samples  is  needed  for  equal  accuracy. 

2.  Analysis  often  required  the  separation  of  data  by 
species,  VSR,  or  by  species  and  VSR  combined.  Since 
these  sample  groups  would  not  be  represented  well  when 
separated  by  these  descriptors,  the  groups  were 
initially  deleted  for  simplicity. 


Table  1.  Distribution  of  sample  trees  by  region  and  by  species. 
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3.  Total  Volume  Estimation 


3 .  1  Introduction 

Tree  volume  is  strongly  related  to  DBH ,  total  height, 
and  form.  Models  developed  for  estimating  total  volume  have 
been  of  three  forms:  (Loetsch  et  al .  ,  1973) 

METHOD  1:  V=  f ( DBH ) 

METHOD  2:  V=  f ( DBH , HT ) 

METHOD  3:  V=  f ( DBH , HT , F ) 


where  V  is  the  volume  from  ground  to  tip 

DBH  is  the  diameter  outside  bark  at  breast  height 
HT  is  the  total  tree  height 
F  is  a  measure  of  form. 

Volume  as  a  function  of  DBH  alone  (Method  1)  is  usually  only 
used  for  small  areas  where  a  reliable  diameter-height 
relationship  can  be  established  and  so,  are  called  local 
volume  equations  (Husch  et  al . ,  1972).  For  larger  areas, 

prediction  of  total  volume  based  only  on  DBH  is  not  as 
accurate  as  the  other  modes  of  prediction;  if  equations  are 
developed  for  a  specific  region  and  species  group,  an 
accurate  function  may  result.  The  more  popular  method  for 
volume  estimation  is  to  estimate  volume  using  height  and  DBH 
(Method  2,  standard  volume  equations).  Heights  must  be 
measured  in  the  sample,  but  this  additional  information 
produces  a  greater  accuracy  of  prediction  (Loetsch  et  al . , 
1973)  . 
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The  third  method  is  potentially  the  most  accurate, 
because  three  parameters  are  included  in  estimating  volume. 
However,  form  measurements  are  difficult  to  obtain  and  may 
be  time  consuming  at  the  field  level  (Burkhart,  1977).  Behre 
(1927)  suggested  that  the  addition  of  form  for  volume 

prediction  could  result  in  functions  useful  for  a  larger 
area,  or  a  greater  number  of  species.  On  the  other  hand, 

form  often  does  not  add  significantly  to  the  estimate 

(Honer,  1965).  Also,  development  of  computer  software  to 
estimate  parameters  using  least  squares  linear  regression 
has  minimized  the  time  and  effort  needed  to  develop  new 
functions  for  localized  areas.  Classification  by  species, 
and/or  area  along  with  Method  2  usually  results  in  a 
sufficiently  accurate  model,  and  so  measurements  of  form  are 
not  necessary.  Because  the  inclusion  of  form  probably  would 
not  significantly  improve  the  volume  estimate  and  the 
measurement  of  form  is  time  consuming,  Method  3  was  not 

tested . 

Many  models  have  been  developed  to  predict  volume  as  a 
function  of  DBH  or  DBH  and  height.  An  assessment  of  existing 
models  is  needed  to  determine  which  model  is  most  accurate 
in  predicting  total  volume  for  the  sample  area. 

Hypothesis  1 

Within  a  method  of  volume  prediction,  no  one  model 
is  consistently  more  accurate  than  other  models. 

Once  models  of  the  two  methods  are  assessed,  the  local  and 

standard  volume  function  methods  should  be  compared  to 
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determine  which  method  is  most  accurate  for  the  sample  data. 


Hypothesis  2 

Models  based  on  DBH  and  height  are  not  significantly 
more  accurate  than  models  based  on  DBH  alone. 

The  way  to  stratify  the  data  to  obtain  the  best  estimate  is 

also  of  interest.  Differences  between  species,  within  a  VSR 

must  be  determined  so  that  categorization  of  data  best 

minimizes  the  variance  of  total  volume. 


Hypothesis  3 

There  is  no  significant  difference  in  total  volume 
between  species  once  variation  due  to  diameter, 
height,  and  regional  differences  is  controlled. 

The  control  of  diameter,  height,  and  regional  differences 

helps  to  reduce  variation  due  to  these  factors,  and  eases 

the  interpretation  of  results  for  species  differences. 

Differences  between  regions  for  a  particular  species  must  be 

determined . 

Hypothesis  4 

There  is  no  significant  difference  in  total  volume 
between  regions  once  variation  due  to  diameter, 
height,  and  species  differences  is  controlled. 

If  Hypotheses  3  and  4  are  not  rejected,  for  practical 

purposes,  one  volume  function  could  be  used  for  all  species 

and  regions. 

All  hypotheses,  together,  are  designed  to  obtain  the 
best  estimate  of  total  volume  using  models  previously 
developed,  and  to  reduce  the  error  of  this  estimate  by 
stratifying  the  total  volume  data. 
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3.2  Methods 

3.2.1  Assessment  of  the  Two  Methods  of  Total  Volume 
Prediction 

Using  the  literature  as  a  guide,  two  models  were  chosen 

to  predict  volume  as  a  function  of  DBH . 

MODEL  1  V  =  a  +  b  DBH  +  c  DBH2 
(Hohenadl  and  Krenn,  in  Loetsch  et  al . ,  1973) 

b 

MODEL  2  V  =  a  DBH 

(Berhout,  in  Loetsch  et  al . ,  1973) 

'transformed  to: 

log  V  =  log  a  +  b  log  DBH  (Husch,  1963). 

where  DBH  is  the  diameter  outside  bark 
at  breast  height(cm) 

H  is  the  total  tree  height(m) 

V  is  the  total  tree  gross  volume(m3) 
a,  b,  c  are  coefficients. 

These  models  were  chosen  because  they  were  representative  of 
the  two  general  types  of  models.  The  first  is  a  simple 
quadratic  model.  The  second  is  a  logarithmic  transformation 
of  variables  used  to  reduce  the  heterogeneity  of  variance 
(Neter  and  Wasserman,  1974).  Other  studies  have  shown  that 
the  variance  of  volume  increases  with  increasing  DBH  or 
total  height  (Cunia,  1964;  Honer,  1965),  so  a  logarithmic 
transformation  may  be  beneficial. 

Five  models  were  chosen  to  predict  volume  from  DBH  and 
he i ght . 

MODEL  3  Schumacher's  logarithmic  volume  (1933) 

b  c 

=  a  DBH  H 


V 
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transformed  to: 

log  V  =  log  a  +  b  log  DBH  +  c  log  H 
MODEL  4  Spurr’s  constant  form-factor  (1952) 

V  =  a  DBH 2  H 

MODEL  5  Spurr's  combined-variable  (1952) 

V  =  a  +  b  DBH  2H 

MODEL  6  Honer  (1965) 

DBH2  =  a  +  b 
V  H 


MODEL  7  Spurr’s  logarithmic  volume  (1952) 

b 

V  =  a ( DBH 2  H ) 

transformed  to: 

log  V  =  log  a  +  b  log ( DBH 2 H ) 

Models  3  and  7  are  logarithmic  transformations, 
originally  developed  so  a  least  squares  linear  regression 
procedure  could  be  used  (Spurr,  1952).  As  with  the 
logarithmic  local  volume  function,  the  logarithmic 
transformation  can  also  aid  in  reducing  the  heterogeneity  of 
variance.  Model  4  uses  DBH  squared  times  height,  which  is 
reportedly  the  best  single  independent  variable  for  the 
prediction  of  volume  (Honer,  1965).  Model  5  is  the 
equivalent  of  Model  4  except  that  the  intercept  is  allowed 
to  vary  from  the  origin.  A  reciprocal  transformation  of 
variables  was  used  by  Honer  (Model  6)  to  obtain  a  more 
homogenous  variance.  The  model  is  relatively  new  and  was 
chosen  because  it  is  associated  with  functions  for 
prediction  of  merchantable  volume.  Models  3,  4,  5,  and  7 
have  been  tested  extensively  with  other  data  sets  and  have 
proved  accurate  (Spurr,  1952;  Schumacher,  1933;  Ellis, 
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1978)  . 

All  seven  models  were  tested  for  accuracy  of  volume 
prediction  using  the  multiple  linear  regression  subprogram 
of  the  Statistical  Package  for  the  Social  Sciences  (SPSS) 
computer  software  package  (Nie  et  a/.,  1975).  Data  were 
separated  by  species  to  reduce  variation  caused  by  species 
differences,  as  well  as  to  introduce  the  possibility  that  a 
model  was  consistently  better  for  one  species  than  for 
another.  Because  the  regression  of  logarithmic  models 
results  in  a  report  of  the  standard  error  of  the  logarithm 
of  volume,  the  standard  error  of  volume  was  approximated 
using : 

SUM  ((actual  -  predicted  volume)2) 

SQRT  - 

error  degrees  of  freedom  (5) 

(Spurr,  1952) 

where : 

SQRT  is  the  square  root 

actual  volumes  are  the  section  tree  volumes 

predicted  volumes  are  obtained  by  taking  the  antilog 
of  the  equation  produced  by  linear  regression 
of  the  transformed  model. 

This  is  only  an  estimate  of  the  standard  error  of  volume, 
but  it  can  be  used  in  a  rough  comparison  to  standard  errors 
from  other  models.  A  similar  procedure  was  used  by  Faurot 
(1977).  The  standard  error  for  Model  6  was  also  approximated 
in  this  manner. 

Models  3  and  7  were  also  fit  using  a  nonlinear 
algorithm  (British  Columbia  Forest  Service,  1981)  for  one  of 
the  species.  The  estimated  coefficients  using  the  nonlinear 
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fit  of  the  multiplicative  model  versus  the  linear  fit  of  the 
logarithmic  transformation  were  compared. 

Hypothesis  1  was  assessed  by  examining  regression 
results  of  models  within  a  particular  method.  The  standard 
error  values,  actual  or  calculated,  were  compared.  Within 
the  logarithmic  models,  the  adjusted  R2  values  were 
compared;  similarly,  a  comparison  of  the  adjusted  R2  values 
for  the  two  linear  models  was  made.3  With  large  sample 
numbers,  the  adjusted  R2  value  approaches  the  R2  value. 

The  logarithmic  version  of  volume  as  a  function  of  DBH 
(Model  2),  and  Schumacher's  logarithmic  volume  function 
(Model  3)  were  compared  as  a  test  of  Hypothesis  2.  An  alpha 
level  of  .05  was  used  to  determine  whether  the  addition  of 
the  logarithm  of  height  significantly  improved  the  model. 

3.2.2  Species  and  Region  Differences 

The  regression  method  for  analysis  of  covariance  (Neter 
and  Wasserman,  1974)  was  used  in  testing  Hypothesis  3. 
Models  were  developed  with  species  and  regions  represented 
as  indicator  variables.  For  species  differences,  the 
following  model  was  formed: 

V  =  a  +  b  DBH 2  H  +cZ,+dZ2+eZ3 

+  f  Z4  +  g  Z5  +  h  DBH  2HZ ,  +  i  DBH  2HZ  2 
+  j  DBH2HZ3  +  k  DBH2HZ4  +  1  DBH  2HZ  5 

3  The  adjusted  R2  value  is  calculated  by  subtracting  the 
ratio  of  the  error  variance  to  the  variance  of  the  total 
from  one  (Fowler  and  Bigelow,  1979). 
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(6) 

where  Z,  =  1  for  white  spruce,  else  0 
Z  2  =  1  for  black  spruce,  else  0 
Z  3  =  1  for  balsam  fir,  else  0 
Z  4  =  1  for  pine,  else  0 
Z  5  =  1  for  aspen,  else  0 
All  Zi  =  0  for  balsam  poplar 

DBH  is  the  diameter  at  breast  height (cm) 

H  is  the  total  tree  height. 

For  differences  in  region,  the  following  model  was  used 


V  =  a  +  b  DBH  2H  +  c  X,  +  d  X2  +  e  X3  + 
f  DBH  2 HX ,  +  g  DBH 2  HX  2  +  h  DBH  2HX3 

(7) 

where  X,  =  1  for  VSR  3,  else  0 

X2  =  1  for  VSR  2,  else  0 

X3  =  1  for  VSR  4,  else  0 

All  Xi  =  0  for  VSR  8. 

Both  equations  were  fit  using  multiple  linear  regression. 
DBH  squared  times  height  was  used  as  a  covariate  to  account 
for  variation  due  to  diameter  and  height  differences.  Hence, 
remaining  variation  was  attributed  to  differences  in  form 
between  VSRs  or  species.  For  differences  in  species,  data 
were  separated  by  VSR  and  vice  versa,  for  easier 
interpretation  of  results. 

Once  the  equations  were  fit,  Scheffe's  method  for 
contrasts  (Neter  and  Wasserman,  1974)  was  used  to  test 
Hypotheses  3  and  4.  The  method  was  applied  to  differences  in 
slope  and  differences  in  intercept  coefficients,  for 
contrasts  between  all  possible  pairs  of  species  and  then 
between  all  possible  VSR  pairs.  For  this  method,  the 
difference  between  coefficients  is  significant  when: 
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B  -  B 

i  j 

- -  >  S 

SQRT  (s2  +  s2  ) 

B  B 

i  j 

(Neter  and  Wasserman , 1 974 ) 

where : 

S2  =  ( k-  1  ) F 

( k- 1 ) , ( n-k- 1 ) , alpha 
SQRT  is  the  square  root 
B  ,  B  are  coefficients  to  be  compared 

i  j 

k  =  the  number  of  treatments 
s 2  s 2 

B  ,  B  are  the  coefficient  variances. 

i  j 


A  confidence  level  of  95%  was  chosen  as  the  test  level;  a 
level  lower  than  this  was  not  considered  accurate,  whereas  a 
level  higher  than  this  would  not  be  practical.4  Although 
Scheffe's  method  is  conservative  relative  to  other  tests,  it 
is  suitable  for  unequal  numbers  in  a  category,  and  many 
contrasts  can  be  tested  (Kleinbaum  and  Kupper,  1978). 


3.3  Results 

3.3.1  Local  and  Standard  Volume  Functions 

The  regression  ,  results  for  local  volume  functions  are 
listed  in  Table  2a,  and  results  for  standard  volume 
functions  are  in  Table  2b.  All  coefficients  were 
significant,  except  the  coefficient  associated  with  DBH  in 
Model  1  for  balsam  fir  and  white  spruce. 

4An  alpha  level  of  .05  was  chosen  for  most  tests  completed 
in  this  project. 
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For  both  of  the  local  volume  functions  tested,  all  R2 
values  were  greater  than  85  percent.  The  approximate 
standard  errors  calculated  for  the  logarithmic  model  (Model 
2)  were  consistently  higher  than  the  actual  standard  errors 
of  the  quadratic  model  (Model  1).  Plots  of  volume  versus  DBH 
indicated  that  variance  of  volume  increased  with  increasing 
DBH,  as  anticipated  (Cunia,1964;  Zar,1968).  The  logarithmic 
transformation  yielded  a  more  homogeneous  variance  of 
volume,  and  therefore  a  more  uniformly  distributed  residual 
plot.  Using  a  stepwise  linear  regression  procedure,  for  the 
quadratic  model,  DBH2  entered  before  DBH  for  all  species; 
the  addition  of  DBH  to  the  regression  was  not  significant 
for  balsam  fir  and  white  spruce. 

Each  of  the  five  models  based  on  volume  as  a  function 
of  DBH  and  height  were  evaluated  using  multiple  linear  least 
squares  regression.  The  comparison  of  standard  errors 
(actual  and  calculated)  indicated  that  Model  7  had  the 
lowest  values  for  each  species.  The  ranking  of  the  remaining 
models  from  lowest  to  highest  standard  errors  was  Model  5, 
then  Model  4,  followed  by  Models  3  and  6,  together.  The 
examination  of  residual  plots  with  Spurr's  two  linear  models 
indicated  a  deparature  from  linear  regression  assumptions, 
in  that  the  variance  of  the  residuals  was  not  homogeneous. 
The  logarithmic  transformation  of  Models  3  and  7  helped  to 
provide  a  more  homogeneous  variance.  Model  6  had 
comparatively  low  adjusted  R2  values,  however  this  value 
applies  to  DBH  squared  over  volume,  rather  than  volume  or 
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the  logarithm  of  volume. 

The  nonlinear  fit  of  Model  3  for  aspen  accounted  for 
94.381  percent  of  the  variance  (R2)  and  for  Model  7,  94.368 
percent  of  the  variance.  The  R2  value  was  assumed  to  be 
comparable  to  the  adjusted  R2  value  reported  for  the  linear 
regressions,  because  of  the  large  sample  number.  The 
standard  error  for  Model  3  was  .10376,  and  for  Model  7, 
.10381.  A  comparison  of  coefficients  follows: 


Model  3 

a  b 

c 

linear 

.0000808  1.97319 

0.72159 

nonlinear 

.00005216  1.97396 

0.87669 

Model  7 

a 

b 

1 inear 

.0000377 

0.98432 

nonlinear 

.00004345 

0.97091 

For  Model  3,  the  coefficients  from  the  linear  fit  of  the 


transformed 

model  were 

close 

in 

value  to  those 

from  the 

nonlinear  fit  of  the 

untransformed  version.  The 

power 

coef  f ic ients 

( ' b T  and 

'c'  )  ’ 

were 

most  similar. 

The  !a? 

coefficient 

was  greater 

with 

the 

linear  fit 

of 

the 

transformed 

model,  indicating 

that 

the  logarithmic 

model 

possibly  predicts  a  higher  volume  than  does  the  nonlinear 
model.  However,  the  graph  of  actual  to  predicted  volumes 
using  the  linear  fit  of  the  model  was  not  visibly  different 
than  the  graph  of  the  actual  to  predicted  volumes  using  the 
nonlinear  version  (Appendix  III).  With  Model  7,  the  power 
coefficients  were  again  similar,  and  the  comparison  of 
graphs  of  the  two  models  showed  no  visible  difference. 
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3.3.2  VSR  and  Species  Differences 

The  results  of  tests  using  Scheffe's  method  for  paired 
contrasts  of  slope  and  intercept  coefficients,  representing 
VSR  and  species  differences,  are  listed  in  Tables  3  and  4. 

Table  3  shows  that  four  of  six  pairs  of  species  in  VSR 
3,  six  of  ten  pairs  in  VSR  2,  six  of  fifteen  pairs  in  VSR  4, 
and  two  of  ten  pairs  in  VSR  8  had  significantly  different 
slope  or  intercept  coefficients.  White  spruce  was  different 
from  pine  and  aspen  in  all  but  VSR  8.  Black  spruce  was 
different  from  pine  and  aspen  in  the  two  Roc ky-Clearwa ter 
regions  (VSR  3  and  2).  Balsam  poplar  was  different  from  pine 
and  aspen  for  the  three  VSRs  tested. 

Table  4  is  a  list  of  S  values  for  differences  in  total 
volume  between  VSR  pairs,  within  a  species  group.  Four  of 
six  pairs  of  VSRs  for  white  spruce,  four  of  six  pairs  for 
black  spruce,  four  of  six  pairs  for  pine,  two  of  six  pairs 
for  aspen,  and  none  of  the  three  pairs  for  balsam  poplar 
were  significantly  different  in  total  volume.  VSRs  2  and  3 
were  not  significantly  different  for  any  of  the  species 
tested.  VSRs  3  (Rocky,  high  elevation)  and  4  (Centre)  were 
different  for  white  and  black  spruce.  VSR  3  and  VSR  8  differ 
for  black  spruce  and  pine.  VSR  2  and  4  were  different  for 
all  species  except  balsam  poplar.  For  pine,  white  spruce, 
and  black  spruce,  the  Rocky  area,  low  elevation  (VSR  2), 
varied  from  the  Lac  La  Biche  area  (VSR  8).  Finally,  VSR  4 
was  different  from  VSR  8  for  white  spruce,  pine,  and  aspen. 


Table  3.  Differences  in  total  volume  between  species  within  a  VSR. 
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3.4  Discussion 

3.4.1  Local  Volume  Functions 

Of  the  two  models  assessed,  the  standard  error  values 
were  lower  for  Model  1;  the  comparison  of  these  standard 
errors,  though,  is  not  unbiased.  The  residual  pattern  for 
Model  1  indicated  a  nonhomogeneous  variance,  a  departure 
from  linear  regression  assumptions.  The  logarithmic 
transformation  helped  to  reduce  the  variance  for  larger 
volumes  and  so  the  resulting  residual  patterns  indicated  a 
more  homogeneous  variance.  The  logarithmic  model  would, 
then,  be  the  recommended  model,  even  though  the  calculated 
standard  error  was  higher  than  the  standard  error  for  Model 
1,  the  quadratic  model.  The  heterogeneity  of  variance  must 
be  ignored  if  Method  1  is  to  be  employed.  A  nonlinear  fit  of 
Model  2  would  be  desirable,  because  the  actual  rather  than 
the  transformed  values  are  minimized  in  the  regression 
procedure . 

3.4.2  Standard  Volume  Functions 

Standard  volume  functions  tested  were  of  three  forms: 

1.  Linear  (Models  4  and  5). 

2.  Logarithmic  or  reciprocal  transformations  to  a  linear 
form  (Models  3,  6  and  7). 

3.  Nonlinear  (Models  3  and  7). 

Comparing  the  standard  error  values,  the  linear  models  have 
been  proven  accurate  in  volume  prediction.  Model  4  which 
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forces  the  regression  through  zero  had  similar  standard 
errors  and  adjusted  R2  values  to  Model  5  in  which  the 
intercept  was  allowed  to  vary.  The  unexpected  higher 
adjusted  R2  values  reported  for  Model  4  were  probably  the 
result  of  the  forcing  the  model  through  zero  (Gordon,  1981). 
The  standard  errors,  however,  showed  the  expected  result  in 
that  Model  5  yielded  slightly  lower  standard  error  values. 
The  advantage  to  the  model  forced  through  zero,  is  that  the 
expected  result  of  zero  volume  with  zero  DBH  and  height  is 
predicted . 

Of  the  transformed  models,  Honer's  model,  which  uses  a 
reciprocal  transformation,  had  low  adjusted  R2  values, 
although  the  adjusted  R2  values  were  reported  for  the  ratio 
of  DBH  squared  over  volume.  The  calculated  standard  errors 
for  Model  6  were  similar  to  those  calculated  for  Model  3. 
However,  Honer's  model  does  require  more  manipulation  of 
variables  to  obtain  the  volume  estimate  than  does  any  of  the 
remaining  models  tested,  and  does  not  provide  a  great  deal 
more  accuracy.  Model  7,  Spurr’s  logarithmic  model, 
consistently  yielded  higher  adjusted  R2  and  lower 
approximate  standard  error  values  than  did  Model  3 
(Schumacher's).  The  model  approximates  a  cone  shape  for  the 
tree  (DBH  squared  times  height),  with  additional  variation 
in  form  being  accounted  for  by  the  variation  in 
coefficients. 

Both  nonlinear  models  produced  high  R2  values  and  low 
standard  error  values  for  aspen,  the  only  species  tested. 
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The  R2  values  were  almost  identical  for  the  two  models.  The 
nonlinear  version  of  Model  3  should  be  slightly  better  in 
estimating  volume  as  an  additional  parameter  is  allowed  to 
vary  than  with  Model  7.  The  proximity  of  results  “for  the  two 
models  then,  indicates  that  very  little  is  gained  by 
allowing  the  variation  of  both  power  coef f ic ients . 5  Further 
testing  of  the  two  nonlinear  models  for  other  species  was 
not  attempted.  However,  the  indication,  from  testing  the 
logarithmic  transformations  of  the  models  is  that  both 
models  yield  accurate  results  over  all  species. 

The  disadvantage  of  linear  models  as  a  group  is  that 
volume  has  unequal  variance,  and  hence  the  assumptions  for 
linear  regression  are  not  met  (Honer,  1965).  The  use  of 
weighted  least  squares  regression  has  been  suggested 
( Cun ia , 1 964 ) ;  the  problem  with  this  is  that  weights  must  be 
known  or  estimated  (Kleinbaum  and  Kupper,  1978).  If  the 
problem  of  unequal  variance  is  overlooked,  Models  4  and  5 
yield  high  R2  values.  Indeed  many  authors  have  used  these 
functions  to  estimate  volume  (Burkhart,  1977;  Myers  and 
Edminster,  1974;  Beagle,  1974;  Bruce  and  Demars,  1974). 

The  use  of  transformations  of  variables  is  advantageous 
in  that  a  more  homogeneous  variance  can  result  (Zar,  1968), 
and  also  a  least  squares  regression  method  can  be  used 

5The  logarithmic  transformations  behave  differently  in  that 
Model  7  provided  more  accurate  results  than  Model  3.  This  is 
because  the  logarithmic  transformation  changes  the  number  of 
dependant  variables  for  Model  3  from  1  to  2  variables.  The 
results  expected  with  the  nonlinear  versions  are  then  not 
applicable  to  the  transformed  versions  of  these  same  models. 
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(Schumacher,  1933).  However,  the  square  of  the  difference 
between  transformed  actual  and  predicted  dependent  variables 
is  minimized  rather  than  the  untransformed  values  (Zar, 
1967).  Also,  the  standard  error  that  is  reported  is  the 
standard  error  of  the  transformed  variable,  and  hence 
comparisons  and  confidence  limits  are  difficult  to  interpret 
(Cunia,  1964).  The  standard  error  of  volume  can  be 
approximated,  however,  some  bias  is  introduced  in  this 
calculation  (Spurr,  1952).  The  empirical  R2  values  could 
also  be  calculated  for  the  transformed  models  and  used  for 
compar i son . 

The  nonlinear  fit  of  Spurr 's  or  Schumacher's  model  is 
the  superior  mode  of  volume  estimation.  For  this  study,  the 
results  for  aspen  are  very  accurate.  Also,  the  advantages  of 
using  the  nonlinear  fit  of  models  are: 

1.  The  nonconstant  variance  of  volume  is  not  a  problem. 

2.  The  square  of  actual  minus  the  predicted  volumes  is 
minimized,  rather  than  the  transformation  of  values. 

3.  The  reported  standard  error  is  the  standard  error  of  the 
volume  (Husch  et  al . ,  1972;  Zar,  1968). 

The  increasing  availability  of  computer  programs  to 
approximate  the  nonlinear  fit  has  made  this  solution  more 
feasible  (Husch  et  al . , 


1972)  . 
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3.4.3  Comparison  of  Local  and  Standard  Volume  Functions 

The  applicability  of  each  of  the  two  methods  of  volume 
estimation  was  assessed.  The  addition  of  height 
(Schumacher's  model)  to  the  logarithmic  model  of  volume  as  a 
function  of  DBH  was  significant  over  all  species.  Hypothesis 
2  was  therefore  rejected,  and  so  the  result  that  volume  as  a 
function  of  DBH  and  height  is  a  significantly  better  method 
of  volume  prediction  was  concluded.  The  addition  of 
polynomial  terms  into  the  regression  of  volume  predicted  by 
DBH,  would  have  accounted  for  a  greater  portion  of  the 
variation,  however  the  returns  dimish  with  each  higher 
polynomial  (Kleinbaum  and  Kupper,  1978).  Even  with  the 
quadratic  model,  for  two  of  the  six  species  tested,  the 
addition  of  DBH,  after  DBH2,  was  not  significant.  Also,  the 
addition  of  height  permits  the  use  of  the  equation  for  areas 
of  differing  height-diameter  relationships.  In  total,  Method 
1  accounts  for  a  large  portion  of  the  variance,  but  for 
large  land  bases,  volume  inventories  based  on  diameter  and 
height  are  suggested  (Husch  et  al . ,  1972;  Loetsch  et  a/., 
1974)  . 


3.4.4  VSR  and  Species  Differences 

Hypothesis  3  could  not  be  wholly  rejected, 
species  pairs  were  different  in  total  volume 
the  remaining  three  VSRs,  approximately  half  of 
pairs  varied  in  total  volume.  Hypothesis  4  was 
white  spruce,  black  spruce,  and  pine.  However, 


as  very  few 
in  VSR  8.  In 
the  species 
rejected  for 
for  aspen  and 
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balsam  poplar,  the  conclusion  that  regions  are  different  in 
total  volume  could  not  be  inferred.  Spurr  (1952)  concluded 
that  species  differences  were  more  significant  than  regional 
differences.  For  this  data  set  it  appears  that  pairs  of 
regions  tested  were  more  significant  in  defining  variation 
in  total  volume.  The  use  of  Scheffe's  method  gives 
conservative  results,  but  the  method  is  commonly  used,  and 
the  magnitude  of  values  produced  indicate  that  similar 
results  would  be  obtained  using  other  methods.  Schumacher 
(1933)  concluded  that  differences  in  total  volume  between 
species  were  mostly  due  to  different  diameter-height 
relationships;  DBH  squared  times  height,  as  a  covariate 
term,  greatly  reduced  the  variation  of  volume. 


3.5  Conclusion 

Overall,  standard  volume  functions  are  more  accurate 
than  local  volume  functions.  Local  volume  functions  could  be 
adequate  for  small  locales  of  homogeneous  diameter-height 
combinations.  Of  the  standard  volume  functions,  the 
nonlinear  models  are  the  first  choice  for  volume  estimation 
(Ellis,  1978;  Husch  et  a\ . ,  1972).  Transformation  of  these 
nonlinear  models  for  use  with  a  linear  regression  method 
also  gives  accurate  results.  If  the  assumption  of  equal 
variance  is  ignored,  linear  models  can  be  used  with  accurate 


results . 
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Once  the  variation  in  volume  due  to  different  diameter 
and  height  combinations  is  accounted  for,  very  little 
variance  remains  and  hence  few  species  and  regions  were 
different  in  total  volume.  Based  on  the  paired  differences, 
the  following  breakdown  should  be  used  to  minimize  the 
variance  of  total  volume  most  efficiently.  Generally,  VSRs  2 
and  3  can  be  combined  for  estimation.  In  all  VSRs,  the  two 
spruce  species  should  be  separated  from  pine.  Aspen  should 
be  estimated  separately  from  other  species;  estimates  made 
for  VSR  4  or  VSR  8  can  be  used  in  the  Rocky  area,  although 
VSR  4  estimates  should  not  be  used  for  VSR  8  and  vice  versa. 
Balsam  poplar  should  be  estimated  separately  from  other 
species  including  aspen,  although  one  volume  function  can  be 
used  for  all  of  the  three  regions 
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4.  Merchantable  Volume 


4 .  1  Introduction 

Merchantable  volume  is  the  portion  of  the  main  stem 
tree  volume  within  specified  utilization  limits.  Separate 
functions  can  be  developed  for  each  utilization  standard, 
although  these  functions  are  often  incompatible  (Burkhart, 
1977;  Belanger  and  Cleroux,  1973).  An  equation  which 
estimates  merchantable  volume  fo*r  many  standards, 
simultaneously,  is  more  desirable.  Functions  for  prediction 
of  merchantable  volume  must  be  flexible  for  a  wide  range  of 
utilization  standards  and  must  predict  total  volume  equal  to 
that  predicted  by  the  total  volume  function  being  used. 

Generally,  two  types  of  simultaneous  estimation  have 
been  used  (Cao  et  al . ,  1980;  Loetsch  et  al . ,  1973).  The 
first  involves  multiplying  total  volume  by  a  factor  less 
than  one  (merchantable  ratio  (MR)).  The  factor  is  estimated 
using  regression  analysis  where: 


Merchantable  Volume 

MR  =  -  =  f  (DBH  ,  H  ,  h  ,  d  ,  .  .  .  ) 

Total  Volume 

where : 

DBH  is  the  diameter  outside  bark(cm)  at  breast 
height (1.3m) 

H  is  the  total  tree  height  (m) 

h  is  the  merchantable  height  (m) 

d  is  the  diameter  inside  bark(cm)  at  'h'. 

The  second  method  is  the  mathematical  integration  of  a  taper 

equation  over  the  merchantable  portion  of  the  tree  (Max  and 
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Burkhart,  1976). 

Because  the  merchantable  ratio  is  multiplied  by  total 
volume  to  obtain  merchantable  volume,  the  estimation  of 
merchantable  volume  using  this  method  is  dependent  on  the 
accuracy  of  both  the  merchantable  ratio  and  total  volume 
function.  Also,  the  regression  estimate  of  the  merchantable 
ratio  rather  than  merchantable  volume  is  obtained.  The 
advantage  of  merchantable  ratio  functions  is  that  they  are 
easily  integrated  into  an  existing  system  for  predicting 
total  volume  (volume-based  system) ,  and  generally  the 
equations  are  not  complex. 

The  advantage  of  the  taper  function  method  is  that  any 
combination  of  lower  and  upper  limits  can  be  used,  and  also 
taper  models  can  be  used  independently  of  total  volume 
functions  (Martin,  1981).  In  addition,  the  merchantable 
height,  merchantable  volume,  total  volume,  and  top  diameter 
are  all  estimated  using  one  model.  The  disadvantage  of  taper 
functions  is  that  they  are  often  complex.  In  addition,  if  a 
volume-based  system  is  already  established,  the  adoption  of 
taper  functions  may  require  the  development  of  a  totally  new 
system  for  volume  estimation  (Cao  et  al . ,  1980). 

An  assessment  of  representive  models  for  each  of  the 
two  methods  to  determine  which  are  appropriate  for  Alberta 
is  needed. 

Hypothesis  6 

Within  a  particular  method,  there  is  no  significant 
difference  in  prediction  accuracy  between  models 
tested . 
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The  most  accurate  of  the  taper  models  and  merchantable  ratio 
models  can  then  be  compared. 

Hypothesis  7 

■There  is  no  significant  difference  in  the  accuracy 
of  the  two  methods  for  merchantable  volume 
estimation . 

The  accuracy  of  the  models  in  predicting  the  merchantable 
volume,  as  well  as  the  computational  ease,  and  compatibility 
with  the  present  system  must  all  be  considered  in  the 
assessment . 

4.2  Methods 

Representative  models  were  needed  for  testing.  Five 
models  were  chosen  which  predict  the  merchantable  ratio. 


1.  Honer's  height  ratio  (1964). 


MR  =  a  +  b  HTOP  +  c  (HTOP) 2 

HT  (HT) 2 


2.  Honer's  squared  diameter  ratio  (1964). 


MR  =  a  +  b  ( DTOP ) 2  +  c  (DTOP) 4 

( DBH ) 2  (DBH) 4 


3.  Honer's  squared  diameter  ratio,  corrected  for  stump 
height  (1974). 

hs  DTOP2  hs  DTOP2  2 

MR  =  a  +  b  ((1+ — )( - — ■))  +  c  ((1+ — )  (- — —  )) 

HT  DBH2  HT 


DBH2 


4. 


Honer's  diameter-height  ratio  (1974). 


4  1 


DTOP  HTOP  2 

MR  =  a  +  b  -  +  c  (  1  -  - ) 

DBH  HT 


5.  Alberta  Forest  Service,  modified  Honer's 
( Lowe ,1981). 


hs  2  DTOP  4  hs  2  DTOP  4  2 

MR  =  a  +  b((1-  —  )-( -  )  )  +  c  (  (  1  —  —  )-(—  )  ) 

HT  DBH  HT  DBH 

where : 

HTOP  is  the  height  above  ground  to  merchantable 
1  imi  t 

HT  is  the  height  from  ground  to  tip 
DBH  is  the  diameter  outside  bark  at  breast  height 
DTOP  is  the  diameter  inside  bark  at  the 
merchantable  limit 
hs  is  the  stump  height. 


Each  model  was  assessed  using  multiple  linear  regression  by 
species.  Almost  all  of  the  functions  tested  were  developed 
by  Honer ;  the  models  also  are  related  in  that  Model  3  is  an 
adaptation  of  Model  2  to  vary  stump  height.  However,  few 
other  models  have  been  developed,  and  Honer's  model  has  been 
tested  with  other  data  sets  (Ker,  1974;  Cao  et  al . ,  1980).  A 
nonlinear  model  developed  by  Burkhart  (1977)  was  shown  by 
Cao  and  others  (1980)  to  be  more  accurate  than  Honer's 
height  ratio  and  diameter-height  ratio  functions  (Models  1 
and  4).  Burkhart's  model,  though,  is  intrinsically  nonlinear 
and  so  a  reliable  nonlinear  regression  algorithm  is 
required.  More  important,  the  model  does  not  have  a 
provision  for  changing  stump  heights;  the  only  independent 
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variable  is  minimum  diameter  at  the  top  divided  by  DBH. 
Because  of  these  reasons,  Burkhart's  model  was  not  tested. 

The  mathematical  integration  of  taper  functions,  the 
second  method,  was  also  tested  using  multiple  regression 
analysis.  Two  models  were  chosen  based  on  accuracy  and 
useability.  The  first  model  chosen  is  a  model  developed  by 
Kozak  and  others  (1969).  The  model  is  basically  a  quadratic 
equation,  conditioned  so  that  when  merchantable  height 
equals  total  height,  the  predicted  squared  diameter  ratio 
(d2/DBH2)  is  zero.  This  model  was  previously  used  by  the 
British  Columbia  Ministry  of  Forests,  and  was  just  recently 
replaced  by  a  segmented  polynomial  equation  (Demaerschalk 
and  Kozak,  1977;  Demaerschalk,  1981;  Kozak  et  al . ,  1969). 


MODEL  6  Kozak  and  others  (1969) 


d 2  h  h 2 

-  =  5(-  -  1 )  +  c ( — —  -  1  ) 

DBH2  H  H2 

h  = ( - bH ) -  SQRT  ( (bH) 2~4c(aH2~d2H2 ) /DBH 2 ) 

2c 

h2-hs2  h3-hs3 

MV  =  kDBH 2  (a(h-hs)  +  b( - )  +  c(— — )  ) 

2H  3H 2 

where  a  =  -b  -c 

a,  b,  c  are  coefficients 

k  is  .00007854  for  estimates  in  metric  units 
hs  is  the  stump  height 
h  is  the  merchantable  height 
H  is  the  total  tree  height, 
d  is  the  minimum  diameter ( i . b . ) 

MV  is  the  merchantable  volume. 

A  model  developed  by  Omerod  (1973)  was  also  tested.  This 
model  is  nonlinear,  conditioned  so  that  zero  is  predicted  as 
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the  minimum  top  diameter  when  the  merchantable  height  equals 
the  total  height.  When  the  merchantable  height  is  breast 
height  (1.3  metres),  the  minimum  top  diameter  predicted  is 
equal  to  the  DBH.  The  model  was  transformed  using  logarithms 
so  that  a  linear  regression  procedure  could  be  employed. 

MODEL  7  Omerod  (1973) 

2b 

H  =(H  -  h ) / ( H  -  1.3) 

DBH2 

transformed  to: 

log  H  =  2b  log ( (H-h)/(H- 1.3)) 

DBH2 

1/b 

h  =  H  -  ( ( d/DBH )  (H  -  1.3)) 

2b+ 1  2b+ 1 

( kDBH 2 ) ( 1 . 3  -  H)  (h  -  H)  (hs  -  H) 

MV  =  -  X  ( -  -  - - - - - -  ) 

( 2b+ 1 )  ( 1 .3  -  H)  ( 1 .3  -  H) 

All  variables  are  as  those  in  Model  6,  above. 

For  both  models,  then,  the  squared  diameter  ratio  is 
predicted  by  a  function  expressing  tree  taper.  This  function 
is  reorganized  to  obtain  merchantable  height,  and  integrated 
for  merchantable  volume  (MV).6 

Multiple  linear  regression  analysis  by  species  was 
implemented  to  determine  coefficients  for  the  functions 
which  estimate  d2/DBH2.  These  coefficients  were  then  used  to 
determine  the  estimated  merchantable  height  and  merchantable 
volume.  The  calculated  merchantable  height  (actual)  was  then 

"Merchantable  heights  must  only  be  measured  to  fit  the 
equation  for  estimating  the  squared  diameter  ratio. 
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regressed  against  the  estimated  merchantable  height 
(predicted)  by  species,  to  determine  how  well  the  fitted 
equation  predicted  merchantable  heights.  Similarily,  the 
merchantable  volume,  actual  and  predicted,  was  regressed  by 
species . 

Testing  of  other  models  describing  tree  taper  could  not 
be  accomplished  within  the  time  period  allowed  for  this 
project.  More  complex  models,  such  as  the  segmented 
polynomial  model  proposed  by  Max  and  Burkhart  (1976)  were 
not  selected  for  testing  because: 

1.  Less  complicated  models,  if  accurate,  are  advantageous. 
If  one  or  '  both  of  the  models  tested  here  proved 
accurate,  the  adoption  of  more  complicated  models  is  not 
practical . 

2.  Segmented  models  require  the  identification  of 

inflection  points,  which  may  be  difficult  to  define. 

Because  of  the  initial  calculations  of  merchantable 
volume,  volumes  used  for  testing  all  models  were  variable 
only  in  top  diameter.  Stump  heights  were  not  varied  (See 
Section  2.2.2).  The  merchantable  volumes  used  in  the  study 
were  then  representative  of  the  present  AFS  stump 
utilization  standard.  Also,  only  those  trees  of  2.5  metres 
in  length  and  greater,  within  the  specified  limits,  were 


included . 
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4.3  Results 


4.3.1  Merchantable 

Ratio 

Method 

Merchantable 

ratio 

models 

as  a  function  of 

variables 

including  merchan 

table 

height 

(Models  1  and 

4)  were 

consistently  more  accurate  (higher  R2 ,  lower  standard  error) 
than  models  which  do  not  include  merchantable  height  as  a 
variable  (Table  5).  Model  3,  which  is  a  modification  of 
Model  2  to  include  stump  height,  did  not  estimate 
merchantable  ratio  substantially  better.  Models  3  and  5  were 
similar  in  accuracy  of  prediction  over  all  species  groups. 

Models  1,  2,  and  4  do  not  include  stump  height  as  a 

variable.  The  fitted  equations,  then,  are  useful  only  in 
predicting  merchantable  volume  for  the  fixed  stump  height 
represented  by  the  original  sample;  a  change  in  stump  height 
is  not  reflected  in  the  predicted  merchantable  ratio.  The 
fitted  equations  for  Models  3  and  5  yielded  different 
merchantable  ratio  values  for  changes  in  stump  height; 
however,  Model  3  did  not  respond  to  changing  stump  heights 
if  the  minimum  top  diameter  was  set  to  zero. 

Each  model  was  tested  to  determine  if  a  value  of  one 
was  predicted  for  total  volume  (merchantable  ratio  of  one). 
Model  1--A  value  of  .92991  was  predicted  as  the  lowest 
merchantable  ratio  for  black  spruce  and  .96354  predicted 
for  pine  as  the  highest  value. 

Model  2--The  range  of  predicted  values  was  .9737  for  white 


spruce  to  1.0048  for  black  spruce. 
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Model  3 --The  range  of  predicted  values  was  .9740  for  white 
spruce  to  1.0045  for  black  spruce. 

Model  4--A  value  of  .9830  was  estimated  for  white  spruce  to 
a  value  of  1.0063  for  pine. 

Model  5--The  range  of  predicted  values  was  .9939  for  white 
spruce  to  1.0235  for  black  spruce. 

The  significance  of  the  variation  from  one  was  not  tested; 
the  fact  that  the  models  did  not  predict  an  exact  value  of 
one  was  considered  only. 

4.3.2  Taper  Models 

Kozak's  model  predicted  merchantable  volume  with 
greater  accuracy  (higher  R2)  than  did  Omerod's  model  (Model 
2)  over  all  species  (Table  6).  However,  Omerod's  model  was 
more  accurate  in  estimation  of  merchantable  height,  and  was 
generally  equal  to  Kozak’s  model  in  prediction  of  the  square 
of  minimum  diameter  (i.b.)  over  DBH .  All  regressions  yielded 
R2  values  of  greater  than  90  percent. 

The  prediction  of  total  volume  using  Kozak's  model 
ranged  from  an  R2  value  of  .96646  and  a  standard  error  of 
.06526  for  balsam  poplar,  to  an  R2  value  of  .98643  and  a 
standard  error  of  .02704  for  balsam  fir.  For  Omerod's  model, 
the  lowest  R2  value  was  .96455  for  balsam  poplar,  with  a 
standard  error  of  .06434;  the  highest  was  .98595  for  balsam 
fir  with  a  standard  error  of  .02751.  Both  models 
overestimated  total  tree  volume  in  large  trees  of  white  and 


black  spruce. 


Table  6.  Results  of  regression  of  taper  models  by  species. 
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4.4  Discussion 


4.4.1  Merchantable  Ratio  Models 

All  models  predicted  the  merchantable  ratio  with  R2 
values  of  over  85  percent.  Honer's  two  models  which  include 
merchantable  height  as  a  predictor  (Models  1  and  4),  gave 
better  results  than  models  without  this  value.  However, 
because  merchantable  height  is  difficult  to  measure  in 
sampling  (Husch  et  a/.,  1972),  these  models  are  not 

practical  unless  an  improved  means  of  estimating 
merchantable  height  is  assumed. 

The  AFS  model  (Model  5)  generally  yielded  high  R2 
values  and  has  the  added  advantages  that: 

1.  Merchantable  height  is  not  required  for  estimation. 

2.  Changes  in  stump  height  are  reflected  by  a  change  in  the 
predicted  merchantable  ratio,  even  when  the  minimum  top 
diameter  is  zero. 

The  disadvantage  to  this  model  was  that  total  volume  is 
generally  overestimated  (predicted  merchantable  ratio 
greater  than  one).  The  reason  for  this  overestimation  could 
be  that  Model  5  is  based  on  a  paraboloid  shape  for  tree 
volume,  whereas  the  original  calculation  of  total  volume 
from  the  sample  volume  employs  a  cone  formula  for  the  top 
section  (See  Section  2.2.1).  The  volume  of  the  top  section 
could  be  overestimated  with  this  model. 

Model  3  yielded  R2  values  similar  to  Model  5.  The 
model,  generally,  underestimated  total  volume  and  did  not 
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respond  to  changing  stump  heights  when  the  minimum  top 
diameter  was  set  to  zero.  The  need  for  a  model  that  predicts 
a  different  merchantable  ratio  for  differing  stump  heights 
with  a  top  diameter  of  zero,  may  not  now  be  necessary, 
although  in  the  future,  as  utilization  standards  change, 
this  feature  may  be  desirable. 

Model  2  gave  results  similar  to  Model  3,  though  stump 
height  was  not  included  as  a  variable  in  the  model.  Any 
change  in  stump  height  then  would  require  the  model  to  be 
fit  for  the  new  data. 

4.4.2  Taper  Models 

The  R2  values  for  all  regressions  were  over  90  percent 
for  both  models  and  all  species.  The  results  indicate,  then, 
that  both  of  the  models  are  accurate.  The  reason  for  these 
high  R2  values  may  be  that  the  variability  of  the 
merchantable  volume  was  limited;  no  stump  heights  were 
varied,  and  only  four  different  minimum  top  diameters  were 
represented.  The  original  testing  of  Model  6  by  Kozak  and 
others  (1969)  showed  that  the  model  accounted  for  95  percent 
of  the  variation  of  the  squared  diameter  ratio.  In  testing 
Model  7,  Omerod  (1973)  suggested  that  the  use  of  a  segmented 
model  yielded  better  results  than  this  simple  model.  Both 
Cao  and  others  (1980)  and  Martin  (1981)  found  that  these  two 
models  were  adequate,  but  not  as  accurate  as  more  complex 
equations.  For  other  estimates,  such  as  the  merchantable 
volume  of  small  log  pieces,  or  the  diameter  below  stump 
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height,  these  models  have  been  found  to  be  biased 
(Demaerschalk  and  Kozak,  1977).  The  two  models  are  accurate, 
then,  in  estimating  the  merchantable  volume,  merchantable 
height,  and  the  squared  diameter  ratio  as  defined  for  this 
study.  For  radically  different  utilization  specifications, 
or  for  small  log  portions,  some  bias  may  result,  although 
further  analysis  is  needed. 


4.4.3  Comparison  of  the  Two  Methods 

A  comparison  of  the  two  methods,  merchantable  ratio 
prediction  and  integration  of  taper  models,  indicates  that 
taper  models  have  the  capacity  for  greater  accuracy  than 
merchantable  ratio  models.  In  terms  of  R2  values,  Hypothesis 
7  can  be  rejected,  and  the  conclusion  that  taper  models  are 
more  accurate  should  be  made. 

Taper  functions  can  be  used  to  replace  a  volume-based 
system;  total  and  merchantable  volume  functions  can  be 
replaced  by  a  single  function  describing  stem  taper.  Even 
simple  models,  such  as  those  tested  in  this  project,  can  be 
more  accurate  than  merchantable  ratio  models  for  predicting 
the  merchantable  volume.  The  taper  functions  are  also 
accurate  in  predicting  total  volume. 

More  complex  taper  functions  will  possibly  be  needed  as 
utilization  standards  change  and  as  the  need  for  precise 
estimates  of  diameter,  merchantable  height,  or  merchantable 
volume  rises  (Bruce  et  al . ,  1968).  Simpler  functions  have 


been  shown  to  be  biased  for  the  base  and  top  portions  of  the 
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tree  stem  ( Demaer schalk  and  Kozak,  1977).  Equations  such  as 
Max  and  Burkhart's  (1976)  segmented  polynomial  function  or 
Omerod's  function  for  a  hyperboloid  form  (1981)  should  be 
tested  for  prediction  of  the  present  as  well  as  possible 
future  merchantable  volume  inventory  needs.  The  problem  with 
these  more  complex  functions  is  that  the  calculations  become 
unwieldy.  Hilt  (1980)  developed  a  computer  program  for 
fitting  the  taper  function  he  developed;  Bruce  and  others 
(1968)  suggest  using  a  merchantable  ratio  for  practical  use 
of  the  function. 

If  a  system  incorporating  merchantable  ratio  functions 
is  already  established,  these  functions  can  be  converted  to 
taper  functions.  The  conversion  procedure  as  described  by 
Clutter  (1980)  involves  the  algebraic  reorganization  of  the 
combined  merchantable  ratio  and  total  volume  function  so 
that  the  minimum  top  diameter  becomes  the  dependent  variable 
(taper  function).  Thus,  instead  of  completely  replacing  the 
established  volume-based  system,  the  advantage  of  an 
accurate  taper  function  can  be  added  by  rewriting  the 
merchantable  ratio  model  to  a  taper  model.  The  accuracy  of 
the  resulting  taper  function  must  be  tested. 

Taper  functions  are  often  incompatible  with  existing 
total  volume  functions  (Munro  and  Demaerschalk ,  1974)  and 
any  inventory  system  which  uses  a  total  volume  function  may 
be  costly  to  replace.  Taper  functions  can  be  mathematically 
conditioned  to  be  compatible  with  total  volume  functions 
(Demaerschalk,  1971  and  1972).  However,  merchantable  ratio 
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functions  are  inherently  more  compatible  as  they  are 
dependent  on  total  volume  functions  for  estimating 
merchantable  volume.  The  adoption  of  a  taper  function 
conditioned  to  predict  total  volume  equal  to  the  existing 
equation  may  prove  more  accurate. 


4.5  Conclusion 

Simple  taper  functions  can  provide  accurate  results  in 
predicting  the  merchantable  volumes  represented  in  this 
project.  They  also  predict  total  volume  accurately.  However, 
if  an  existing  system  incorporates  a  total  volume  function 
with  a  merchantable  ratio  model,  the  continued  use  of 
merchantable  ratio  models  is  more  practical  and  estimates 
are  quite  accurate. 

For  future  inventory  analysis  systems,  the  following 
options  are  available: 

1.  A  taper-based  system  could  be  adopted  for  predicting 
total  and  merchantable  volume,  as  well  as  merchantable 
height  and  minimum  top  diameter. 

2.  Existing  merchantable  ratio  models  could  be  converted  to 
taper  models.  The  resulting  taper  function,  though,  may 
not  be  as  accurate  as  other  taper  models  available. 

3.  The  use  of  taper  functions,  conditioned  to  predict  total 
volume  equal  to  the  existing  total  volume  function  may 
be  the  more  accurate  and  practical  route. 

Overall,  merchantable  ratio  functions  and  simple  taper 
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functions  are  accurate  for  prediction  of  the  present  needs 
of  merchantable  volume  estimation.  A  taper  model  is 
advantageous,  but  may  not  be  feasible  for  use  with  an 
existing  inventory  system.  Further  testing  of  both  methods 
to  determine  limitations  is  required. 


5.  Defect  Estimation 


5 . 1  Introduction 

Internal  defect  is  the  amount  of  missing  or  crumbling 
wood.  Each  type  of  wood  product  has  associated  losses 
because  of  internal  defect.  Defect  volumes  between  trees 
vary  greatly  depending  on  the  disease  susceptibility  of  the 
tree,  the  availability  of  pathogens,  and  various  site 
characteristics.  For  this  reason,  estimation  of  defect  is 
difficult;  resulting  estimates  are  often  statistically  poor. 

The  aim  in  estimation  is  to  obtain  the  best  estimate  of 
defect  by  maximizing  the  reduction  in  the  variance  of  the 
estimate,  while  minimizing  the  number  of  discrete  and 
continuous  predictors  employed.  A  constraint  on  these 
predictors  is  that  they  must  be  quickly  and  easily  measured 
on  standing  trees.  Species,  site  indicators,  and  cull 
suspect  indicators  are  examples  of  discrete  variables 
previously  used  (Smith,  1973).  These  variables  can  be 
represented  either  through  preliminary  segregation  of  data 
or  by  indicator  variables  in  a  regression  model.  Continuous 
variables  such  as  DBH,  age,  and  number  of  conks  can  be 
directly  used  as  independent  variables  in  a  regression  model 
(Aho,  1974);  they  may  also  be  translated  to  discrete  units 
by  dividing  the  range  of  the  variable  into  classes. 

For  the  AFS  tree  section  data,  defect  volumes  for 
sawlog  and  pulp  product  types  were  calculated  (See  Section 
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2.2.3).  VSRs,  species,  and  cull  suspect  classes  are  the 
discrete  variables  which  can  be  used  as  predictors  of 
defect.  To  best  use  these  variables  to  reduce  the  variance 
of  defect,  statistically  significant  differences  must  be 
determined.  First,  differences  in  defect  volume  between 
species  and  between  VSRs  must  be  defined. 

Hypothesis  8 

There  is  no  significant  difference  in  percent  defect 
(pulp  or  saw)  between  species. 

Hypothesis  9 

There  is  no  significant  difference  in  percent  defect 
(pulp  or  saw)  between  VSRs. 

Significantly  different  groups  should  be  separated  before 
estimating  defect  values.  They  may  alternatively  be 
represented  as  indicator  variables  in  a  regression  model. 

The  following  four  defect  indicator  classes  (cull 
suspect  classes)  are  presently  used  by  the  AFS: 

1  .  Nonsuspect 

2.  Scars  and  Others 

3.  Old  Broken  Top 

4.  Conk  and  Punk  Knot. 

The  importance  of  these  classes  as  indicators  of  internal 
defect  was  of  interest. 

Hypothesis  10 

There  is  no  significant  difference  in  percent  defect 
(pulp  or  saw)  between  cull  suspect  classes. 

Again,  significantly  different  cull  classes  can  be 
represented  as  indicator  variables  in  a  regression  model. 
The  expected  result  is  that  the  'nonsuspect’  class  would 
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indicate  less  defect  than  the  other  classes. 

The  continuous  variables  represented  in  this  study  were 
age  and  DBH.  The  relationship  between  age  and  defect,  and 
DBH  and  defect  was  determined  by  testing  the  following 
hypotheses : 

Hypothes i s  1 1 

There  is  no  significant  relationship  between  age  and 

percent  def ect (pulp/sawlog ) . 

Hypothesis  12 

There  is  no  significant  relationship  between  DBH  and 

percent  def ect (pulp/sawlog ) . 

If  DBH  or  age  is  proven  significant  in  estimating  defect, 
data  can  be  separated  into  classes,  and  an  average  defect 
percent  calculated,  or  these  variables  can  be  represented 
directly  as  independent  variables  for  regression  estimation. 

Once  the  significant  discrete  and  continuous  predictors 
are  found,  a  regression  model  to  incorporate  all  of  the 
pertinent  predictors  can  be  developed. 

%defect  =  f (DBH , age , cull  indicators,.  .  .  ) 

If  all  significant  discrete  variables  are  represented  as 
indicator  variables  in  the  regression  model,  the  model  could 
become  large.  Hence,  some  of  the  discrete  variables  should 
be  used  to  categorize  data  before  regression  analysis  is 
attempted.  Even  though  all  significantly  different  variables 
are  represented  either  in  the  regression  model  or  in  the 
segregation  of  data,  the  resulting  regression  model  may  be 
statistically  insignificant. 
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Hypothesis  13 

A  regression  model  to  estimate  percent  defect  will 
not  be  significant. 

Percent  defect  is  the  dependent  variable  in  all  tests  rather 
than  actual  defect  volume.  By  using  percent  defect,  the 
variance  in  defect  volumes  contributed  by  differences  in  the 
total  tree  volumes  is  controlled. 


5.2  Methods 

Hypotheses  8  through  10  were  tested  using  subprogram 
ANOVA  of  the  SPSS  statistical  package,  with  an  alpha  of  .05. 
Differences  between  pairs  of  group  means  were  assessed  using 
Scheffe’s  method  for  multiple  comparisons.  ANOVA  were 
completed  for  the  following: 

1.  Percent  defect  with  species,  by  VSR  (Hypothesis  8). 

2.  Percent  defect  with  VSR,  by  species  (Hypothesis  9). 

3.  Percent  defect  with  cull  class,  by  VSR  species,  and  by 
VSR/species  group  (Hypothesis  10). 

Two  way  analysis  of  variance  to  simultaneously  test 
Hypotheses  8  and  9  was  not  attempted,  because  of  the 
limitations  of  the  SPSS  version  used;  paired  contrasts  could 
not  be  done  with  two-way  analysis  of  variance  (Nie  et  al . , 
1975).  To  replace  this  two-way  analysis  of  variance,  data 
were  separated  by  species,  when  testing  for  VSR  differences, 
and  by  VSR,  when  testing  for  species  differences.  Similarly, 
data  were  separated  by  species,  by  VSR,  and  by  VSR/species 
groups  in  testing  Hypothesis  10. 
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Regression  of  percent  defect  with 
separately,  was  employed  to  determine  if 
variables  significantly  improve  the  est 
percent  (Hypotheses  11  and  12).  An  alpha 
used  to  determine  whether  the  resulting 
significant.  Data  were  separated  by  VSR  and 
The  results  of  all  tests  were  then  in 
regression  model  was  developed.  The  reg 
tested  (Hypothesis  13)  using  multiple  1 
analysis,  with  a  significance  level  of  .05. 


age  and  DBH , 
these  continuous 
imate  of  defect 
level  of  .05  was 
regression  was 
species. 

terpreted  and  a 
ression  model  was 
inear  regression 


5.3  Results 

5.3.1  Species  Differences 

Differences  in  percent  pulp  defect  between  species  were 
significant  in  all  but  VSR  3  (Table  7).  Aspen  differed  in 
percent  pulp  defect  from  all  other  species  in  VSRs  2  and  4. 
In  VSR  8,  aspen  differed  from  black  spruce  and  pine.  No 
other  species  pairs  were  significantly  different. 

For  all  VSRs  tested,  differences  in  percent  saw  defect 
between  species,  as  a  group,  were  significant.  Defect  in 
aspen  differed  from  all  species  in  VSRs  2  and  4,  and  from 
all  species  except  balsam  poplar  in  VSR  8.  Balsam  poplar 
defect  was  different  from  pine  in  VSRs  2,  4,  and  8,  and  from 
black  spruce  in  VSRs  2  and  8. 


Table  7.  Differences  in  percent  def ect ( pu 1  p/saw )  between  species  by  VSR. 
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In  VSR  8 

( Lac 

La  Biche- 

Athabasca 

region ) ,  balsam 

poplar 

di  f  f  ered 

in 

percent  saw 

defect  fi 

rom  white  spruce. 

In  VSR  3 

( Roc  ky , 

high 

elevation ) 

,  black 

spruce  and  pine 

defect 

percents 

were 

different . 

5.3.2  VSR  Differences 

Differences  in  percent  pulp  defect  between  VSRs  as  a 
group  were  significant  for  two  of  five  species  tested  (Table 
8).  VSR  4  differed  from  other  regions  for  aspen.  The  VSRs, 
as  a  group  differed  in  pulp  defect  for  balsam  poplar, 
although  the  difference  is  too  small  for  any  significant 
differences  between  individual  regions. 

In  three  of  five  species,  differences  in  percent  saw 
defect  between  VSRs  were  significant.  For  defect  in  black 
spruce,  VSR  3  (Rocky,  high  elevation)  differed  from  other 
VSRs.  Pine  in  Rocky,  low  elevation  (VSR  2)  differed  from 
pine  in  Lac  La  Biche-Athabasca  region  (VSR  8).  As  with  pulp 
defect,  VSR  4  differed  in  saw  defect  for  aspen  from  other 
VSRs. 

5.3.3  Defect  Indicator  Class  Differences 

Cull  class  differences  by  VSR  for  pulp  defect  (Table  9) 
were  significant  in  all  VSRs  tested.  Class  4  differed  from 
other  classes  in  VSRs  2,  4  and  8,  and  class  1  differed  from 
2  in  VSRs  3  and  8.  For  saw  defect,  differences  were  noted 
between  cull  class  4  and  other  classes  in  VSRs  2  ,  4  and  8. 


Table  8.  Differences  In  percent  def ect ( pu 1  p/saw )  between  VSRs  by  species. 


62 


V) 

0) 

> 

> 

u 

f— 

p— 

1-  c 

■H 

E 

c 

■e 

as 

c 

O 

0 

C 

U  L 

(0 

£_ 

E 

r 

0 

0) 

0 

4- 

0 

*F 

o 

lii  4- 

**- 

CO 

L. 

4- 

4- 

•F 

4- 

■F  CO 

4- 

u.  -  c 

■*- 

C  13 

c 

•p. 

•D  0) 

c 

0  C 

w 

0  T3 

c 

ai  d) 

O) 

L  0 

L. 

a  C 

OS 

■w  3 

•f- 

0 

0 

0  0 

■«- 

Q  C  +> 

in 

4-  - 

4- 

4- 

in 

0  d)  in 

-H 

4-  ^r 

4- 

4-  CO 

■e 

o  n  t- 

in  c 

•r- 

-r- 

•»- 

in  c 

£-  0 

T3  - 

T3 

■0  - 

£-  0 

4-  (0 

•*-  C- 

CN 

CN 

t_ 

a 

0  0 

CO 

CN  00 

0  0 

s  c 

a  4- 

in 

in 

a  4- 

0)  a. 

4- 

a  oi 

Ql  Q1 

01  01 

4- 

<  ~  in 

o  — 

CO  CO 

00  l/) 

00  00 

O  — 

in  > 

Z  *0 

>  > 

>  > 

>  > 

z  ts 

in 

m  o 

CO 

* 

* 

* 

o 

a  — 

O’ 

G) 

CD 

CO 

r- 

in  -h 

T— 

h* 

in 

ID 

>  in 

T- 

CN 

CD 

o 

*- 

—  -H 

IT) 

CO 

—  (0 

■O' 

<  -H 

in 

u. 

in 

l-  <U 

> 

> 

> 

> 

a 

*— 

t — 

o  c 

■F 

■F 

E 

•F 

d» 

c 

c 

c 

0 

c 

LU  L. 

(0 

ns 

0 

L. 

0 

d) 

a 

O 

o 

4- 

o 

Li.  4- 

**■■ 

-r- 

•r— 

4- 

4- 

4- 

4- 

-F  CO 

4- 

UJ  —  C 

f- 

c 

■r~ 

T3  a) 

C 

c 

c 

0  "0 

c 

a  a) 

D) 

OS 

O) 

L.  C 

D) 

■h  3 

-f- 

•*- 

0  0 

•*- 

C 

in 

in 

in 

4- 

in 

ns  (U  in 

+■> 

+■> 

•F 

4-  CO 

0  XI  u 

in  c 

in  c 

in  c 

•t- 

in  c 

t-  d) 

c.  d) 

L.  0 

T3  * 

a  0 

CL  <*.  ns 

•>-  i_ 

•<-  i_ 

L. 

CN 

•*-  L. 

a 

ns  d) 

ns  0 

0  0 

0  0 

-1  c 

a  4- 

a  4- 

0.4- 

in 

a  4- 

as  a 

4- 

4- 

4- 

01  01 

4- 

D  -  in 

0  — 

0  — 

0  — 

(/)  l/) 

o  — 

in  > 

2  TS 

Z  T3 

Z  *0 

>  > 

Z  *0 

CL 

in  o 

cn 

00 

CD 

* 

* 

ct  — 

in 

US 

0) 

^r 

in  +4 

os 

Lf) 

CO 

r- 

>  in 

CN 

O 

6 

o 

T“ 

in 

CO 

—  ns 

CO 

<  +» 

in 

U. 

C- 

<D 

0 

0 

o 

u 

3 

3 

a 

L. 

L. 

0 

a 

a 

CL 

in 

in 

W 

as 

E 

ds 

JL 

c 

0 

o 

V 

U 

0 

0 

in 

as 

■r- 

0 

c 

a 

' — 

a 

r 

»— 

•r- 

in 

0 

in 

3 

CO 

a. 

< 

CO 

0 

> 

0 


in 

O 


d) 

.c 


0 


c 

0 

u 


c 

OS 

in 

* 


Table  9.  Differences  in  percent  def ect ( pu 1  p/saw )  between  cull  classes  by  VSR  group. 
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Classes  1  and  2  were  also  different.  In  VSR  8,  all  pairs  of 
classes  were  different  in  saw  defect. 

Table  10  shows  cull  class  differences  in  percent 
pulp/saw  defect  by  species.  For  pulp  defect,  no  cull  class 
differences  were  significant  for  black  spruce;  all  pairs 
were  different  for  white  spruce.  For  pine  and  aspen,  'conk 
and  punk'  (class  4)  was  different  in  defect  from  other 
classes.  Classes  1  and  2  were  also  different  for  aspen.  For 
balsam  poplar,  cull  class  3  was  different  from  cull  classes 
1  and  2.  For  saw  defect,  cull  class  4  differed  from  other 
classes  for  white  spruce,  balsam  poplar,  pine  and  aspen. 
Classes  1  and  2  varied  in  percent  saw' defect  for  all  species 
except  balsam  poplar.  For  white  spruce  and  balsam  poplar, 
classes  2  and  3  were  significantly  different. 

Table  1 1  is  a  summary  of  ANOVA  tests  for  differences  in 
pulp  defect  between  cull  classes  by  VSR/species  group.  In  11 
of  20  VSR/species  groups,  cull  class  differences  were 
significant.  Generally,  in  aspen,  'conk  and  punk'  (class  4) 
differed  from  'nonsuspect'  (class  1)  and  'scars  and  others' 
(class  2).  For  VSR  2/black  spruce,  and  VSR  8/aspen,  classes 
1  and  2  were  different.  Cull  class  4  differed  from  class  1 
in  VSR  2/balsam  poplar.  For  VSR  4/white  spruce,  class  4 
differed  from  other  classes.  'Scars  and  others ' (class  2)  was 
different  from  other  classes  for  VSR  8/white  spruce.  'Conk 
and  punk  knot'  differed  from  other  classes  for  VSR  8/pine. 
Class  3  was  different  from  classes  1  and  2  for  VSR  8/balsam 
poplar . 


Table  10.  Differences  in  percent  def ect ( pu 1  p/saw )  between  cull  classes  by  species  group. 
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Table  11.  Cull  class  differences  in  percent  pulp  defect  by  VSR/species  group. 
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ANOVA  results  for  cull  class  differences  in  percent  saw 
defect  by  VSR/species  group  are  listed  in  Table  12. 
Differences  were  significant  in  13  VSR/species  groups.  No 
differences  were  significant  for  all  VSR  3/species  groups. 
For  aspen  in  the  remaining  VSRs ,  'conk  and  punk  knot'  (class 
4)  differed  from  other  classes.  Generally,  for  white  and 
black  spruce,  classes  1  and  2  were  different.  Classes  4  and 
1  were  different  from  other  classes  for  VSR  8/  pine.  In  VSR 
4,  class  1  and  3  were  different  for  balsam  poplar,  however, 
few  observations  were  obtained  for  class  3,  'old  broken 
top'  . 


5.3.4  Relationship  of  Age  and  Percent  Defect 

For  13  of  20  VSR/species  groups  tested,  the 
relationship  between  percent  pulp  defect  and  age  (Table  13) 
was  significant  for  an  alpha  of  .05.  The  range  of  R 2  values 
for  significant  regressions  was  .0095  to  .2960;  only  three 
R2  values  were  above  .15.  The  higher  R2  values  resulted  with 
aspen.  For  black  spruce,  the  regression  was  not  significant 
in  any  VSR  tested.  In  three  VSRs,  the  relationship  was  so 
insignificant  that  no  regression  line  could  be  calculated  at 
the  default  tolerance  level  for  SPSS. 
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Because  the  defect  values  were  very  small  or  zero,  a  regression  line  could  n 
calculated  using  the  default  tolerance  levels  of  the  SPSS  regression  package. 
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Regressions  of  percent  saw  defect  with  age  showed 
in  12  of  20  VSR/species  groups,  the  relationship 
significant.  Aspen,  again,  yielded  some  of  the  highes 
values  with  .3052  for  VSR  4.  The  range  of  R2  value 
significant  regressions  was  0.0165  to  .3052;  only  thre 
the  R2  values  were  above  .15. 

5.3.5  Relationship  of  DBH  and  Percent  Defect 

Table  14  is  a  list  of  the  regression  results  for 
and  percent  pulp  and  saw  defect. 

For  13  of  20  VSR/species  groups,  the  regression 
percent  pulp  defect  was  significant.  As  with  age 
regression  appears  to  be  most  significant  for  aspen, 
range  of  R2  values  for  significant  regressions  was  .01 
.2895,  with  only  one  of  the  R2  values  above  .15. 

For  11  of  20  VSR/species  groups,  the  .regression  o 
with  percent  saw  defect  is  significant.  The  range  o 
values  for  significant  regressions  was  .0195  to  .3415 
regression  for  aspen  yielded  a  substantially  higher  R2. 
two  values  were  above  .15. 
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5.3.6  Multiple  Regression  of  Percent  Defect 

Using  results  from  previous  studies  as  a  guide,  and  the 
results  obtained  in  testing  Hypotheses  8  through  12,  the 
following  regression  model  was  developed  for  testing. 

%  defect  =  a  +  bDBH  +  c  AGE  +  d  Z,  +  eZ2 
+  fZ3 

(7) 

where  DBH  is  the  diameter ( o . b . )  at  breast  height(cm) 
AGE  is  the  age  measured  at  stump  height 


1 , 

if 

'nonsuspect',  0  else 

1 , 

if 

'scars  and  others', 

0  el 

1 , 

if 

' old  broken  top' ,  0 

else 

all  Zi  are  zero  for  'conk  and  punk  knot'. 

The  model  is  an  adaptation  of  a  model  used  by  Aho  (1974)  and 
is  similar  to  that  used  by  Brown  (1934).  Coefficients  were 
fitted  using  stepwise  linear  regression  by  species  (Table 
15).  For  pulp  defect,  the  regression  was  significant  for  all 
species  but  black  spruce.  The  regression  for  saw  defect  was 
significant  for  all  species.  The  R2  values  for  both  pulp  and 
sawlog  defect  were  less  than  40  percent,  with  higher  values 
for  aspen.  In  most  cases,  age  entered  first  into  the 
regression;  the  addition  of  DBH  after  age  was  usually 
insignificant  and  sometimes  reduced  the  overall  significance 
of  the  regression. 


Table  15.  Multiple  regression  results  for  predicting  percent  def ect ( pu 1  p/saw ) . 
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5.4  Discussion 


5.4.1  Species  Differences 

For  both  pulp  and  saw  defect,  Hypothesis  8  was 
rejected;  •  species  differences  were  significant.  Pairwise 
comparisons  of  species  revealed  that  aspen  contributed  most 
to  the  species  differences;  aspen  differed  from  other 
species  in  both  pulp  and  sawlog  defect.  Morawski  (1967) 
reported  that  aspen  differed  in  defect  volume  from  other 
species  in  Ontario.  Basham  (1958)  concluded  that  the 
susceptibility  of  aspen  to  decay  hinders  the  use  of  the 
species.  The  results  of  this  study  support  these  earlier 
reports . 

In  addition  to  aspen,  balsam  poplar  differed  from  other 
species  for  sawl'og  defect.  This  indicates  that  deciduous 
species  should  be  estimated  separately  from  coniferous 
species  in  defect  estimation.  As  an  isolated  case,  in  the 
Rocky-Clearwater ,  high  elevation,  region,  pine  and  black 
spruce  differed  in  saw  defect  percentages.  The  reason  for 
this  difference  is  not  obvious.  Possibly,  a  pathogen  may  be 
present  in  this  area  which  affects  one  of  the  species. 
However,  many  other  factors  including  the  hardiness  of  the 
genotypes,  and  the  effects  of  the  environment  on  the  species 
may  be  contributing  to  this  difference  in  defect. 
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5.4.2  VSR  Differences 

Differences  in  pulp  defect  between  VSRs  were  generally 
not  significant  and  so  Hypothesis  9  could  not  be  rejected. 
VSR  4  differed  from  other  VSRs  for  aspen  and  hence  this  VSR 
should  be  segregated  from  other  VSRs  in  estimating  aspen 
defect. 

For  sawlog  defect,  only  half  of  the  six  species  tested 
were  different  in  defect  between  VSRs,  and  so,  again, 
Hypothesis  9  could  not  be  rejected.  For  black  spruce,  VSR  3 
should  be  estimated  separately  from  other  VSRs.  VSRs  2  and  8 
are  different  in  pine  defect;  VSR  4  should  be  separated  from 
other  VSRs  for  aspen. 

VSR  differences  are  mostly  insignificant  except  for  a 
few  isolated  cases.  The  black  spruce  trees  in  the 
Rocky-Clearwater  high  elevation,  were  different  from  black 
spruce  in  other  areas;  many  factors  including  environmental 
stress,  a  difference  in  genotype,  or  the  presence  of 
pathogens  could  be  contributing  to  this  difference.  The 
difference  in  sawlog  defect  for  pine  in  VSR  2  from  VSR  8  is 
probably  explained  by  the  fact  that  two  species,  lodgepole 
pine  ( Pinus  contorta  var.  latifolia)  and  jack  pine  ( Pinus 
banKsiana)  were  represented  (Hosie,  1969).  In  VSR  2, 
lodgepole  pine  is  most  prevalent,  whereas  in  VSR  8,  jack 
pine  is  the  prevailing  pine  species.  Between  these  two  areas 
(for  example  VSR  4),  the  species  hybridize,  so  the 
difference  in  defect  was  not  as  notable  between  VSR  4  and 
VSR  2,  or  between  VSR  4  and  VSR  8.  VSR  4  differs  from  all 


76 


other  VSRs  for  aspen,  indicating  that  the  site  and/or 
species  differ  in  this  region.  Wall  (1971)  indicated  that 
differences  in  defect  between  clones  is  apparent  in  aspen; 
Fowells  (1965)  cited  both  studies  which  supported  and 
studies  which  refuted  the  idea  that  site  is  related  to 
defect  in  aspen.  For  Alberta,  this  study  indicates  that 
differences  in  VSRs  for  aspen  defect  volume  are  significant 
for  some  regions. 

Most  forest  management  agencies  separate  lands  into 
zones  for  estimating  defect.  The  British  Columbia  Ministry 
of  Forests  separate  the  land  base  into  ecological  zones 
(Dobie  and  Kasper,  1974;  British  Columbia  Forest  Service, 
1976).  Manitoba  stratifies  by  management  unit  (Manitoba 
Forest  Inventory,  1981)  and  Saskatchewan  expects  to 
differentiate  by  forest  administration  units  (Lindenas, 
1981).  This  study  indicates  that  separation  by  VSRs  in 
Alberta  does  not  significantly  improve  the  estimate  of 
defect  percent. 

5.4.3  Defect  Indicator  Class  Differences 

In  all  of  the  VSRs  tested,  for  percent  pulp  defect, 
cull  class  differences  were  significant;  3  of  4  VSRs 
differed  in  percent  saw  defect.  Hypothesis  10  was  then 
rejected.  Classes  1  and  4,  and  1  and  2  were  almost  always 
different;  for  a  few  VSRs,  'old  broken  top'  (class  3)  was 
also  different.  By  species,  cull  classes  differed  in  pulp 
defect  for  five  of  six  species  and  in  sawlog  defect  for  all 
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species.  Over  half  of  the  VSR/species  groups  showed 
significant  differences  in  values  for  pulp  and  saw  defect 
between  cull  classes.  Again,  classes  1  and  4,  and  1  and  2 
were  most  commonly  different. 

In  general,  cull  classes  are  significantly  different  by 
species,  by  VSR,  and  by  VSR/species  group.  These  differences 
are  hard  to  quantify;  a  further  complication  is  that  stems 
often  possess  more  than  one  indicator  (Allemdag  and  Honer, 
1972).  The  classes  tested  in  this  study  may  not  be  the  best 
indicators  of  defect,  although  the  use  of  'conk',  'scars', 
and  'knots'  as  indicators  is  common  (Bailey  and  Dobie,  1977; 
Aho ,  1974).  Other  indicators  of  defect  that  have  been  used 

are  'rot  diameter  at  one  foot'  (Brown,  1934),  growth  rate 
(Stayton  et  al  . ,  1970),  and  'frequency  of  unsound  knot' 

(Allemdag  and  Honer,  1972).  However,  the  use  of  'conk'  for 
aspen,  and  'scars  and  others'  for  balsam  poplar  have  been 
found  to  be  the  most  commonly  recorded  indicators  of 
internal  defect  in  Slave  Lake,  Alberta  (Bailey  and  Dobie, 
1977).  The  results  of  this  study  show  that  the  AFS  cull 
‘indicator  classes  significantly  contribute  to  the  estimation 
of  defect. 

5.4.4  Relationship  of  Age  with  Percent  Defect 

Age,  as  an  independent  variable,  significantly  reduced 
the  variance  of  pulp  and  saw  defect  for  over  half  of  the 
VSR/species  groups  tested,  however  the  associated  R2  values 
were  very  low.  Hypothesis  11  was  therefore 


rejected . 


_ 
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Although  Morawski  (1967)  stated  that  cull  was  related  to  age 
for  all  tree  species,  this  study  indicated  that  for  Alberta, 
age  and  defect  are  significantly  related  only  for  a  few 
species,  particularly  aspen  and  pine,  and  these 
relationships  are  relatively  weak.  Bailey  and  Dobie  (1977) 
found  that  age  is  significantly  related  to  defect  for  aspen; 
age  was  not  significantly  related  to  defect  for  balsam 
poplar.  Basham  (1958)  also  concluded  that  defect  increased 
with  age  in  aspen.  The  R2  values  were  higher  for  aspen,  but 
the  age  is  difficult  to  measure  in  this  species;  annual 
rings  are  often  indescr iminant  (Maini  and  Coupland,  1964). 
Age  as  a  predictor  of  defect  is  significant  for  some 
species,  but  may  not  be  easily  measured  and  resulting  R2 
values  are  low. 

5.4.5  Relationship  of  DBH  with  Percent  Defect 

Because  age  is  difficult  to  measure,  DBH  is  often 
substituted  as  a  predictor  (Morawski,  1967).  For  these  data, 
similar  results  were  obtained  for  the  regression  of  DBH  with 
defect  as  with  age  and  defect,  and  hence  Hypothesis  12  was 
also  rejected.  The  range  of  R2  values  was  again  low,  with 
aspen  and  pine  yielding  the  highest  values.  The  use  of  DBH, 
instead  of  age,  is  feasible  though  neither  DBH  nor  age 
greatly  reduces  the  variance  of  defect.  DBH  as  a  predictor 
of  defect  for  aspen  and  pine  appears  to  be  more  significant 
than  in  other  species.  The  reduction  in  variance  with  other 
species  is  minimal. 


79 


5.4.6  Multiple  Regression  of  Percent  Defect 

Multiple  regression  to  estimate  defect  has  been 
attempted  by  many  authors  (Aho,  1974;  Aho  and  Simonski, 
1975;  Stayton  et  al . ,  1970).  The  results  obtained  through 
multiple  regression  of  the  model  which  uses  DBH,  age  and  AFS 
cull  indicator  classes  had  a  maximum  R2  value  of  less  than 
40  percent  for  pulp  and  saw.  In  the  stepwise  regression 
procedure,  DBH  was  always  entered  after  age.  The  later 
addition  of  DBH  sometimes  reduced  the  significance  of  the 
overall  regression,  and  the  variable  DBH  was  not  significant 
given  the  other  variables  were  in  the  model.  Other  variables 
such  as  ’rot  diameter  at  one  foot  above  ground’, 
successfully  used  by  Brown  (1934),  might  be  included  to 
improve  the  regression.  Some  preliminary  stratification  by 
region  may  have  improved  the  regression  estimate,  although, 
differences  in  defect  between  VSRs  were  found  to  be 
significant  only  for  a  few  VSR  pairs. 


5.5  Conclusions 

Because  the  multiple  regression  model  developed  yields 
low  R2  values,  the  defect  percentage  estimates  produced  are 
not  very  accurate.  For  large  scale  inventories,  on  the 
average,  the  estimates  may  be  adequate,  but  for  small  areas, 
and  few  species,  they  would  not  be  accurate. 

If  the  multiple  regression  model  is  to  be  used,  the 
following  breakdown  to  reduce  variance  is  suggested,  based 
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on  the  regression,  and  analysis  of  variance  tests. 

1.  Aspen  must  be  separated  from  other  species  for  both  pulp 
and  sawlog  defect  estimates.  VSR  4  should  be  separated 
from  other  VSRs  for  aspen. 

2.  Balsam  poplar  should  be  also  be  separated  from  other 
species,  though  one  estimate  can  be  used  for  all  three 
VSRs  tested. 

3.  For  VSR  3,  sawlog  defect,  pine  and  black  spruce  should 
be  separated. 

4.  Cull  classes  can  be  employed  to  reduce  defect  variance. 
The  'old  broken  top'  class  rarely  helped  to  reduce 
variance;  another  indicator  class  may  be  more  related  to 
defect . 

5.  Classification  by  age  or  DBH  class  is  suggested  for  some 
species,  particularly  aspen  and  pine. 

The  regression  model  for  estimation  of  defect  volumes  using 
discrete  and  continuous  predictors  of  defect  is  significant, 
but  hardly  accurate.  Other  indicators  of  defect  for  use  in 
the  regression  may  improve  the  estimate.  Another  method  of 
estimating  defect  such  as  a  mill  study,  or  local  cull  study 
may  be  required  for  a  more  accurate  estimate  of  defect 
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6.  Final  Conclusions  and  Recommendations 

Methods  for  estimating  total,  merchantable  and  defect 
volumes  for  four  regions  of  Alberta  were  assessed.  Although 
many  methods  have  been  developed  for  each  type  of  volume, 
the  methods  tested  were  the  most  popular  because  of  accuracy 
and  useability. 

For  total  volume  estimation,  standard  volume  functions, 
which  use  DBH  and  height  as  predictors,  are  more  accurate 
than  local  volume  functions  (DBH  only).  Of  the  standard 
volume  functions,  the  nonlinear  models  are  the  preferred 
models.  Logarithmic  transformations  of  these  models  to  a 
linear  form  are  also  accurate,  however  transformed  values 
are  then  minimized  in  the  least  squares  procedure.  Linear 
models,  especially  Spurr's  constant  or  variable  form-factor 
equations,  are  accurate,  but  the  assumption  of  equal 
variance  is  not  met.  VSR  differences  appear  to  be  more 
significant  than  species  differences,  once  variation  due  to 
diameter  and  height  is  accounted  for. 

Merchantable  volume,  as  specified  in  this  project,  is 
better  predicted  through  the  integration  of  taper  models, 
than  through  the  use  of  merchantable  ratio  models.  Both 
models  produce  very  adequate  results.  The  choice  of  which  of 
the  two  methods  to  implement  depends  on: 

1.  The  present  inventory  system.  If  the  system  is  rigidly 
based  on  a  total  volume  function,  taper  models  may  be 
conditioned  to  be  compatible,  although  merchantable 
ratio  models  may  be  more  practical  for  this  situation. 
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2.  The  type  and  precision  of  estimates  which  will  be 
required  for  future  inventories.  If  prediction  of 
merchantable  height,  merchantable  volume,  and  minimum 
top  diameter  is  required,  for  varied  portions  of  the 
stem,  taper  models  are  more  suited. 

The  choice  of  methods  for  merchantable  volume  estimation,  as 
defined  in  this  study,  remains  a  decision  of  what  method  is 
most  practical  for  use. 

Defect  volumes  must  be  estimated  separately  by  product 
type.  The  wide  variability  of  defect  makes  the  development 
of  an  adequate  regression  model  difficult,  and  the 
regression  model  developed  was  of  marginal  value.  Other 
variables,  more  correlated  with  defect,  may  be  found,  or 
another  type  of  defect  estimation  such  as  a  mill  study  could 
be  employed.  The  regression  model  developed  can  yield  a 
general  estimate  of  defect  percentages  for  use  in  a  large 
scale  inventory.  To  obtain  the  best  estimate  of  defect  using 
the  regression  method  and  present  variables,  data  should  be 
separated  by  AFS  cull  indicator  classes.  Age  or  DBH 
significantly  reduce  the  variability  of  some  species, 
particularly  aspen.  Aspen  should  be  estimated  separately 
from  other  species,  as  should  balsam  poplar.  Some  regions 
can  be  combined,  for  example,  one  estimate  of  percent  defect 
(pulp  or  saw)  for  balsam  poplar  can  be  used  in  the  three 
VSRs  tested. 

Other  species  and  VSRs  in  Alberta  must  be  assessed  for 
each  of  the  three  types  of  volume  analysis.  A  nonlinear 
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regression  model  for  total  volume  prediction  is  most 
accurate  if  a  time  and  cost  efficient  statistical  package 
for  nonlinear  regression  is  available.  Additional  taper 
models  should  be  examined  for  use  in  Alberta.  For  defect, 
the  regression  method  for  large  scale  inventories  is 
promising,  however,  other  variables  for  predicting  this 


value  are  needed 
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APPENDIX 


APPENDIX 


I.  Sample  tally  sheets  showing  parameters  recorded 
in  tree  sectioning. 


IA. 


IB. 


Record  sheet  in  use  prior  to 


Record  sheet  for  1977  through 


977  . 


to  1982. 
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APPENDIX  |  A 

Alberta  Forest  Service  —  Inventory  Survey  Section 
TREE  SECTION  TALLY  SHEET 


(Use  reverse  for  diagrams  and  comment)  Transator 
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92 


<j\ 


r— 1 

CO 

1 

m 

CN 

cr> 

rH 

« 

• 

CO 

CN 

CO 

m 

o 

o 

Du 

Du 

cu 

r-H 

o 

2 

CO 

2 

2 

N 

\ 

Du 

E-t 

r— I 

Eh 

y 

CO 

>H 

-} 

H-l 

< 

E-* 

O 

3 

z 

o 


CO 

w 

C4 

a; 

E-* 


CO 

Cl, 

< 


w 

£h 

< 

Q 


W 

5 

Z 


- 


APPENDIX  II 


Flowchart  of  Computer  Programs 
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APPENDIX  II 


Flowchart  of  Computer  Programs 


APPENDIX  III.  Graphical  comparison  of  nonlinear  and 
linearized  versions  of  Models  3  and  7  for  estimating  tot 


volume . 


APPENDIX  I IIA.  Model  3. 
APPENDIX  I  I  IB.  Model  7. 
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APPENDIX  IV.  Selected  statistics  for  individual  volume 

estimation  in  Alberta. 


Table  IV. 1.  Regression  coefficients  for  VSR/species  groups  using  the 
logarithmic  transformation  of  Schumacher's  volume  function. 
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Aspen  Ave.  %defect  19.88  13.24  27.27 

No.  Trees  22  2  85 

Balsam  poplar  Ave.  %defect  13.31  17.67  30.93 


*■ 


$ 


